NIC  FILE  COPY  AD  All  8744 


AFGL-TR-82-0020 


TECHNIQUE  DEVELOPMENT  FOR  WEATHER  FORECASTING 


Alan  M.  Gerlach  (ed.) 


Systems  and  Applied  Sciences  Corporation 
6811  Kenilworth  Avenue 
P.  O.  Box  308 
Riverdale,  MD  20737 


30  November  1981 


Final  Report  for  Period  1  December  1980  -  30  November  1981 


Approved  for  public  release;  distribution  unlimited 


AIR  FORCE  GEOPHYSICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
HANSCCM  AFB,  MASSACHUSETTS  01731 


& 


i' 


82  u  8  3  0  08 1 


Qualified  requestors  may  obtain  additional  copies  from  the 
Defense  Technical  Information  Center.  All  others  should 
apply  to  the  National  Technical  Information  Service. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (W> ten  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


1  REPORT  NUMBER 

AFGL-TR -8 2-0020 


mWM 


4.  TITLE  (mid  Submit) 


TECHNIQUE  DEVELOPMENT  FOR  WEATHER  FORECASTING 


author^; 


Alan  M .  Gerlach  (ed.) 


PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Systems  and  Applied  Sciences  Corporation 
6811  Kenilworth  Avenue,  P.  O.  Box  308 
Riverdale,  MD  20737 


11  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Final  Report 
1  Dec  80  -  30  Nov  81 


6  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERfaJ 


F 1962 8-81 -C-00 39 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


61102F,  62101F 
6670C  0AB 


12  REPORT  DATE 


Air  Force  Geophysics  Laboratory 
Hanscom  AFB,  MA  01731 

Charles  Burger /Manager/LY _ 


14.  MONITORING  AGENCY  NAME  ft  AODRESSi'/J  dillerent  from  Controlling  Oilier)  IS.  SECURITY  CLASS,  (ol  this  report) 


30  November  1981 


13.  NUMBER  OF  PAGES 

187 


UNCLASSIFIED 


I5».  DECLASSIFICATION-  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (ol  the  abstract  entered  in  Block  20,  il  dlllerent  from  Report) 


19.  KEY  WORDS  (Continue  on  reverse  aide  it  neceaaary  and  identity  by  block  number) 

_  *  * 

GLOBAL  MOISTURE  DISTRIBUTION  SATELLITE  IMAGERY  ANALYSIS 

INITIALIZATION  CLUSTERING 

MESOSCALE  FORECASTING  CLOUD  BRIGHTNESS 

INTERACTIVE  FORECASTING  (McIDAS)  DOPPLER  RADAR  WIND  RETRIEVAL 


20.  ABSTRACT  (Continue  on  reverie  aide  If  neceaaary  and  identity  by  block  number) 

This  report  summarizes  weather  research  conducted  by  SASC  in  several 
technical  areas:  in  numerical  weather  prediction,  investigation  of  the  phys¬ 
ical  and  spectral  characteristics  of  global  moisture  distribution,  and  exam¬ 
ination  of  the  generation  of  harmonics  for  analysis  of  meteorological  data  in 
the  vertical  and  horizontal;  development  of  McIDAS-based  interactive  aids  for 
formulating  mesoscale  forecasts;  development  of  specialized  computer  programs 
for  sensor  and  technique  evaluation;  in  satellite  meteorology,  interactive 


do  ,:sr»  1473  COITION  OF  1  NOV  6ft  1$  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  Data  Entered) 


TINPT.Agg  TF  TED 


SECURITY  CL ASSlPlC  ATlOH  OP  Tu'*‘  P *GE(Whmn  Dm tm  Entmrmd) 


■*’  x*>> 


FOREWORD 


This  is  the  Final  Report  under  Contract  F19628-81-C-0039  with  the 
Meteorology  Division,  Air  Force  Geophysics  Laboratory.  The  contract 
extended  from  December  1,  1980  to  November  30,  1981. 

The  report  consists  of  narrative  descriptions  of  work  accomplished 
during  the  contract  period  and  two  technical  papers  reporting  research 
achievements.  Work  which  was  reported  on  previously  (in  Scientific  Re¬ 
port  No.  1)  is  referenced  in  the  appropriate  place  but  is  not  repeated. 

The  contract  called  for  effort  in  five  discrete  areas,  reported 
upon  in  the  monthly  R&D  Status  Reports  under  these  headings: 

Numerical  Weather  Prediction 
Mesoscale  Forecasting 
Satellite  Meteorology 
McIDAS  Operation  and  Enhancement 
Radar  Meteorology 

This  Final  Report  is  organized  similarly. 

Personnel  associated  with  these  areas  at  any  time  during  the  contract 

were : 

Numerical  Weather  Prediction  -  Donald  C.  Norquist,  M.S. 

Isidore  M.  Halberstam,  Ph.D. 
Shu-Lin  Tung,  M.S. 

Chris  Johnson,  B.S. 

Joan -Marie  Freni,  B. A. 

Mesoscale  Forecasting  -  Joan  M.  Ward,  A.B. 

Maureen  F.  Hampton,  B.A. 
Russell  C.  Dengel,  B.S. 

Gail  A.  Dengel,  M.s. 

Randy  Schechter,  M.S. 

Charles  F.  Ivaldi,  Jr.,  B.S. 
Michael  E.  Niedzielski,  B.S. 
Scott  D.  Hamilton 

Satellite  Meteorology  -  Gail  A.  Dengel,  M.S. 

Robert  P.  d'Entremont,  B.S. 
Michael  E.  Niedzielski,  B.S. 


McIDAS  Operation  and  -  John  M.  Powers,  B. A. 

Enhancement  _  . 

Barry  A.  Mareiro 

Gail  A.  Dengel ,  M.S. 

Radar  Meteorology  -  F.  Ian  Harris,  Ph.D. 

Chen-Hung  Chang,  M.S. 

Glenn  R.  Smythe,  M.S. 

Principal  Investigator  was  Alan  M.  Gerlach,  Ph.D. 

Subcontracts  with  Harris  Corporation  Computer  Systems  Division  and 
Perkin-Elmer  Customer  Service  Division  provided  on-site  maintenance  ser¬ 
vices  for  the  McIDAS  and  Weather  Radar  Data  Processing  System  respectively. 

Reports  were  prepared  by  the  scientists,  engineers,  and  mathemati¬ 
cians  identified  in  the  Table  of  Contents  with  their  technical  reports. 
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SECTION  1 . 


NUMERICAL  WEATHER  PREDICTION 


A.  A  STUDY  OF  THE  PHYSICAL  AND  SPECTRAL  CHARACTERISTICS  OF  ULOBA1 


MO I STl  1  RE  D I STR I  BUT I ON 

I.  Introduction 

Forecast i ng  atmospheric  moisture  distributions  and  ret;  ions  and  inten¬ 
sities  of  precipitation  remains  one  of  the  leading  problems  in  numerical 
weather  prediction.  While  parameterizations  of  large-scale  and  cumulus- 
scalt  precipitation  within  the  models  have  received  much  of  the  blame  for 
forecast  i nadequacies ,  a  factor  which  cannot  be  overlooked  but  which  lias 
received  too  little  attention  is  the  preparation  of  the  initial  moisture 
:  ieh!s. 

r  .  >'  the  current  moisture  analysis  schemes  used  to  support  o per¬ 

il  t 1 2 ,  i !  models  are  simply  univariate  offshoots  of  the  analysis  procedures 
for  the  other  mass  field  parameters.  Little  attention  is  given  to  the 
similarities  or  differences  that  may  exist  between  moisture  and  other 
mass  field  distributions.  That  the  moisture  analysis  is  as  important  to 
the  forecast  as  the  model  itself  has  been  shown  by  Phillips  (1978)*  among 
others.  He  interchanged  moisture  analyses  in  the  same  model  with  substan¬ 
tially  different  precipitation  forecast  results. 

Another  potential  area  for  problems  is  the  use  of  relative  humidity 
instead  of  absolute  humidity  (the  variable  carried  in  most  models)  as  the 
analyzed  quantity.  In  this  study  we  have  taken  a  first  step  toward  ad¬ 
dressing  the  particular  problem  of  moisture  analyses.  Taking  both  a  sum¬ 
mer  and  a  winter  synoptic  case,  we  have  studied  the  global  distribution 

of  moisture  from  two  perspectives.  First,  using  observed  data  from  the 

2 

First  GARP  Global  Experiment  (FGGE)  of  the  Global  Atmospheric  Research 
Proqram  (GARP),  global  plots  were  subjectively  analyzed  to  facilitate  an 
examination  of  the  fields  in  physical  space.  Second,  using  Hough  analyses 
from  FGGE  as  first-quess  fields  and  performing  successive  corrections  on 
them  using  the  observations,  representations  of  the  observed  fields  at 
grid  points  were  constructed  for  spectral  analyses  of  the  fields.  This 
report  compares  the  results  of  the  mass  and  motion  fields  with  the  moisture 

1.  Phillips,  N.  A.,  1978:  A  test  of  finer  resolution.  Office  Note  171, 

National  Meteorological  Center,  National  Weather  Service,  NOAA. 

2.  Obtained  from  Department  of  the  Air  Force,  OL-A,  USAF  Environmental 

Technical  Applications  Center  (MAC) ,  Federal  Building,  Asheville, 

NC  28801. 
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field,  identifying  similarities  and  differences,  to  determine  the  need 
for  a  separate  analysis  scheme  for  moisture. 

II.  Construction  and  Comparison  of  Observed  Fields 

Two  synoptic  times  were  chosen  for  the  extraction  of  data  from  FGGF 
data  tapes:  January  15,  1978  at  OOGMT  (referred  to  as  JAN) ,  and  July  20, 
1978  at  OOGMT  (JUL) .  These  dates  were  chosen  only  because  they  allowed 
the  greatest  ease  in  extracting  data  from  available  tapes.  Synoptic 
cases  were  chosen  rather  than  a  time-average  of  a  longer  period  of  obser¬ 
vations  because  of  the  objective  to  improve  analyses  for  instantaneous 
forecast  models,  not  climate  prediction  models.  However,  it  is  obviously 
desirable  to  identify  typical  similarities  and  differences  between  the 
fields  for  an  analysis  scheme  usable  for  all  synoptic  situations.  "he 
preferred  method  would  be  to  analyze  many  synoptic  events  separat  and 

then  averaqe  the  statistics,  especially  in  the  spectral  analyse...  ,ee 

this  is  only  a  preliminary  study,  it  is  not  claimed  that  the  rel.‘ 
ships  discovered  are  totally  representative  of  typical  relations, 
tween  the  various  meteorological  fields.  However,  many  of  the  relation¬ 
ships  discovered  in  this  study  would  be  expected  to  hold  up  in  such  an 
expanded  study. 

In  extracting  the  data  from  the  FGGE  tapes  a  cutoff  time  of  ±  two 
hours  was  used.  Only  upper-air  (radiosonde)  data  were  used  for  this 
study.  In  addition,  only  observations  which  satisfied  the  horizontal 
and  vertical  quality  control  checks  were  used.  Once  each  data  set  was 
extracted,  individual  observations  which  were  clearly  erroneous  were  re¬ 
moved  —  an  average  of  about  two  observations  for  each  parameter  per  man¬ 
datory  pressure  level.  The  resulting  data  values  were  plotted  on  larqe 
mercator  projection  maps  and  were  hand-analyzed.  Examples  of  these  maps 
are  shown  in  Figs.  la-1.  Such  maps  were  prepared  for  both  JUL  and  JAN  cases 
and  for  the  parameters  T  (temperature,  1  K)  and  q  (specific  humidity,  g/kg) 
for  all  mandatory  levels  between  the  surface  and  300  mb,  for  Z  (geopoten¬ 
tial  height,  gpm)  at  500  mb,  and  mean  sea  level  (MSI.)  pressure.  While  all 
or  the  resulting  maps  were  studied,  further  discussions  in  this  report 
will  concentrate  on  the  850  mb  and  500  mb  levels. 

Because  of  the  relative  sparsity  of  data  in  the  Southern  Hemisphere 
and  over  all  oceans  except  possibly  the  North  Atlantic,  no  strong  argu¬ 
ments  can  be  made  for  the  accuracy  of  the  contours  drawn  in  these  regions. 
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Only  rawinsonde  and  surface  land/marino  observations  were  used  and  only 
made  within  two  hours  of  OOGMT  (before  and  after)  were  included.  Contoi 
the  oceans  and  Southern  Hemisphere  are  not  necessarily  accurate. 


(Contours  were  not  drawn  to  fit  data  values  which  appeared  suspect  when 
compared  with  surroundinq  values.)  Therefore,  our  discussions  of  physical 
fields  will  focus  on  the  Northern  Hemisphere  and  primarily  over  and  near 
the  continents.  Because  of  the  low  spatial  resolution  inherent  in  such 
a  representation,  only  comparisons  in  the  average  and  large-scale  eddy 
domains  can  be  made.  No  attempt  is  made  to  identify  small-scale  similar¬ 
ities  or  differences  in  these  physical  space  representations. 

Climatologically ,  the  zonally  averaged  moisture  is  a  maximum  at  the 
equator  and  decreases  monotonically  toward  the  poles  (Oort  and  Rasmusson, 
1971).  2  This  is  reflected  in  the  cases  represented  here.  Of  greater 
interest,  however,  are  the  departures  from  the  zonal  mean  which  may  re¬ 
sult  from  (1)  land  masses  and  (2)  disturbances  in  the  zonal  mean  flow. 

The  former  can  be  accounted  for  by  climatological  studies  of  moisture 
distribution.  A  product  of  one  such  study  is  shown,  in  Fig.  2,  where 

4 

Peixoto  (1970)  plots  the  time-averaged  vertically  integrated  precipit- 
ab)e  wjtet  iv  dot  im  >i  m  t  >  i  ms  of  the  time-averaged  specific  humidity  q  as 


-  1  /-300  mb 


•'  r  d  / 

*4f 


where  q  is  the  acceleration  due  to  gravity  and  p  is  pressure.  By  remov¬ 
ing  the  zonal  mean  at  each  latitude  from  both  the  time-averaged  field 
and  any  synoptic  field,  one  can  compare  the  magnitude  of  the  departure 
from  the  zonal  mean  due  entirely  to  the  land  masses  with  that  due  to  ed¬ 
dies  which  appear  in  the  chosen  synoptic  case.  For  example,  when  the 
JAN  q  field  was  used  for  the  synoptic  field,  the  two  sources  of  departures 
from  the  zonal  mean  were  found  to  be  of  the  same  order  of  magnitude.  Cal¬ 
culations  using  q  data  for  a  station  at  43.72  N  and  65.25  W  (near  the 
moisture  maximum  associated  with  the  northwest  Atlantic  storm)  indicated 
a  departure  of  2.3  x  10  kg  m”2  in  precipitable  water  from  the  zonal  aver¬ 
age  for  JAN  at  45°  N  computed  from  Pressman  corrected  q  fields  (explained 
in  Section  III).  This  largest  eddy  departure  is  compared  to  the  largest 

3.  Oort,  A.  H. ,  and  H.  M.  Rasmusson,  1971:  Atmospheric  Circulation 

Statistics.  NO A A  Prof.  Pap.  No.  5  (NTIS-72-50295) .  323  pp. 

4.  Peixoto,  J.  P.,  1970:  Water  vapor  balance  of  the  atmosphere  from  five 

years  of  hemispheric  data.  Nordic  Hydrol . ,  2,  120-138. 
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departure  from  the  zonal  mean  of  the  time-averaged  values  on  Fig.  2  (over 
the  Sahara)  of  -1.8  x  10  kg  m  ^  (where  the  time-averaged  zonal  mean  was 
computed  from  Oort  and  Rasmusson  for  30°  N,  yearly  average). 

Considering  both  contributions,  we  see  perturbations  in  the  zonally 
averaqed  moisture  field  for  both  dates.  Identifying  cyclonic  systems 
from  the  MSL  pressure  map,  we  see  disturbances  over  the  west  and  east 
coasts  of  the  U.S.,  eastern  Europe,  and  off  the  east  coast  of  Asia  for 
JAN.  In  all  four  cases  the  moisture  contours  on  the  850  mb  diagram  are 


Fig.  inl.ii  ster eographic  map  ot  vertically  integrated  time  mean  pro- 

cipitable  water  W  fojj  the  period  May  1,  1^58  -^April  30,  1963.  Units 
for  W  are  x  10  kg  m  (from  Rosen  et  a).,  I  v^7  y  ~l ;  taken  originally  from 
Peixoto,  1970) . 

5.  Rosen,  R.  D. ,  t).  A.  Salstein,  and  J.  P.  Peixoto,  1979:  Variability 
in  the  annual  fields  of  large-scale  atmospheric  water  vapor  trans¬ 
port.  Mon.  Wca.  Rev. ,  107,  26-37. 
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pushed  northward  ahead  of  the  surface  low  and  southward  behind  it.  This 
is  also  true  of  the  major  circulations  which  appear  on  the  JUL  MSL  pres¬ 
sure  map.  The  "Bermuda  Hiqh"  centered  in  the  Atlantic  is  pushing  moist¬ 
ure  into  the  central  U.S.,  while  the  somewhat  weaker  cyclonic  circula¬ 
tion  over  the  Persian  Gulf  area  is  advecting  moist  air  into  central  China 
and  dry  air  into  the  eastern  Mediterranean.  This  correlation  between  flow 
perturbations  and  moisture  perturbations  is  also  evident  at  500  mb,  al¬ 
though  less  so  because  of  the  smaller  moisture  values  at  this  level. 

The  temperature  advection  patterns  are  at  least  as  clearly  seen  as 
those  of  moisture.  Positive  departures  from  the  zonal  mean  exist  where 
southerly  flow  occurs  and  negative  departures  exist  in  regions  of  norther¬ 
ly  flow.  This  is  the  nature  of  the  eddies,  as  they  act  to  offset  the 
pole-equator  temperature  gradient  which  is  primarily  responsible  for  their 
origin  and  maintenance.  Thus  the  large  scale  eddies  bring  cold,  dry  air 
equatorward  and  warm,  moist  air  poleward.  The  correlation  between  warm 
and  moist  areas  and  dry  and  cold  areas  in  the  midlatitudes  is  clear  in 
the  figures. 

An  explanation  should  be  given  as  to  why  specific  humidity  was  used 
in  the  study  rather  than  relative  humidity.  Initially,  maps  such  as  those 
in  Fig.  1  were  prepared  for  relative  humidity  in  the  same  manner  as  for 
the  other  parameters.  It  was  seen  immediately  that  relative  humidity 
observations  showed  a  high  degree  of  spatial  irregularity  —  pockets  of 
hiqh  humidity  were  often  located  next  to  pockets  of  low  humidity  at  all 
levels.  Such  a  discontinuous  field  is  very  difficult  to  analyze  subjec¬ 
tively,  and  the  ability  of  any  objective  analysis  scheme  to  replicate 
such  a  highly  irregular  field  is  questionable.  The  specific  humidity  was 
somewhat  more  irregular  than  other  fields  but  was  more  continuous  spatially 
than  relative  humidity  and  therefore  was  used  in  this  study. 

III.  Construction  and  Comparison  of  Analyzed  Fields 

To  compare  more  closely  the  scales  of  variation  in  the  moisture  fields 
with  those  of  the  other  parameters,  it  was  necessary  to  obtain  representa¬ 
tions  of  the  observed  values  on  a  regular  qrid.  This  could  be  done  through 
a  direct  interpolation  of  the  observations  onto  grid  points,  but  any  choice 
of  interpolation  would  act  to  modify  the  spectral  character  of  the  field. 
Because  this  study  concentrates  on  characteristics  of  observed  fields,  it 
is  important  to  construct  an  analysis  field  which  is  most  representative 


of  the  observations.  The  Flattery  analyses  in  the  Level  III  FGGE  data 
set"  include  the  biases  of  this  particular  analysis  scheme  (e.g.,  a  trun¬ 
cation  at  15°  longitude  wavelength)  and  cannot  be  said  to  be  the  most 
representative  of  the  data.  This  fact  was  iemonstruted  by  Chu  and 
Parrish  (1977)^  when  they  showed  that  a  simple  scheme  that  transfers  data 
to  grid  points  resulted  in  a  humidity  field  that  better  represented  ob¬ 
servations  than  the  Flattery  analysis  (described  by  Cooley,  19746 7 8)  for 
a  case  over  the  eastern  F.S.  The  method  chosen  for  this  study  used  the 
Flattery  analyses  as  first-guess  fields,  then  observations  were  used  to 
correct  the  field  in  successive  scans  following  directly  (in  a  univariate 
manner)  the  well-known  scheme  developed  by  Cressman  (1959)  .  The  sate 

number  of  scans,  scan  i adii,  and  distance  weighting  function  given  by 
Cressman  were  used  in  this  study. 

Upon  completion  of  the  Cressman  correction  fields,  both  the  original 
field  (Ftd)  and  the  collected  field  (CFG)  were  interpolated  bilinearly  to 
the  observation  locations  for  .  comparison  with  the  observed  values  (OB). 
The  following  root  mean  square  error  quantities  were  then  calculated  for 


each  field: 


CFG 


"  N 

’  N  ’ 

.2  (CFG.  -  OB.)2 

s 

2  (FG,  -  0B.)“ 

1-1  i  1 

'  rfg 

1=1  l  l 

N 

N 

where  N  is  the  number  of  observations  used  in  each  field  for  correction 

of  the  first  guess.  The  quantity  1  —  ( R  /R  )  can  then  bo  used  as  a 

(  FG  r  G 

measure  of  how  much  impact  the  observations  had  in  correcting  the  first 
guess  field;  the  closer  to  unity,  the  more  the  corrected  field  is  like 
the  observed  field.  Values  for  this  quantity  are  given  in  Table  1. 

After  their  construction,  the  corrected  fields  were  plotted  on  a 
mercator  projection  map  in  order  to  compare  their  gross  features  with 
the  subiective  analysis  maps.  The  figures  for  q’  and  v1,  departures  from 


6.  Chu,  R. ,  and  D.  Parrish,  1977;  Humidity  analyses  for  operational  pre¬ 

diction  models  at  the  National  Meteorological  Center.  Office  Note 
140,  National  Meteorological  Center,  NWS,  NOAA. 

7.  Cooley,  D.  S.,  1974:  A  description  of  the  Flattery  global  analysis 

method  -  No.  1.  Tech.  Proc.  Bui.  No.  J()5,  NWS,  NOAA. 

8.  Cressman,  G.,  1959:  An  operational  objective  analysis  system.  Mon. 

Woa.  Rev.,  87,  367-174. 


25 


<  >!  moisture  ;  ,er- 


zonal  averages,  are  included 

in  Figs.  !a-d.  The 

1  l  i  < 

!  nun  st  u  j  e  ;  ,er- 

turbation 

mentioned  earlier 

are  seen,  and  a  rti 

lut  i<  >r. 

wit  ii  v  1  is  cl  ear 

tor  those 

areas  from  the  maps. 

TAB LK  1. 

FORMULATION  OF  CKLSFMAN-CORKKCTbb  FI:  Li 

1/15. 

7H 

7  a  o 

Quantity 

No .  of  Obs . 

1  -  (  K  .  k , .  , ) 

i , :  i  ibs 

1  -  (  K  /k 

CFO  p 

t  sre 

0070 

.22!  u  ■ 

.23 : 

. 29298 

T  !  000 

354 

.  ,470t> 

296 

.574  6 .3 

1  . 

601 

.  28  59  1 

•  1  9 

. 36135 

1  7<j<j 

559 

.40271 

*t!~  4 

.  38666 

1  ‘>00 

1,21 

.  205115 

(.44 

.25932 

1  400 

50  K 

. 35490 

r:  3 

. 36533 

1  loo 

57  1 

.  39409 

(24 

.  39354 

>  lo.O 

148 

. 50520 

290 

.59354 

,  6  5  ; 

5  84 

.  457  o8 

0:  )6 

. 50569 

'Ml) 

50< 

.  38845 

526 

.  30987 

.  500 

54  1 

.48644 

554 

.46568 

,  400 

■ .  ,■> 

-t  •  - 

.1,0817 

4.54 

.41448 

VOi' 

124 

.  7142.3 

426 

.41211 

1  MSI. 

2079 

.  22058 

225H 

.23681 

7.  50  0 

f>  15 

.  24800 

,  39 

. 32142 

U  850 

764 

.  .3 1  467 

oil 

.  .32220 

U  500 

678 

. 2104  1 

742 

. 1 2747 

V  850 

764 

.  3 3704 

8  31 

.  34187 

V  500 

678 

.  2507  3 

74  2 

.24275 

The  corrected  fields  were  analyzed  spectrally  using 

both  zonal  and 

spherical 

harmonic  expansion 

;.  The  departure  from  ‘he  zonal  average  A’ 

was  calculated  for  each  2.5° 

latitude  and  a  fast 

Fourier 

transform  (FFT) 

algorithm 

was  used  to  obtain 

the  coefficients  a  , 
m 

b  of 
m 

the  expansion  of 

each  ruant 

i  t-  y  A : 

M 

A1  (^,\)  -  XI  (a  (v)  cos  m\  t  b 

_ ,  V  m  m 

m=  1 

(4)  sir, 

12) 

where  if  is 

latitude,  \  longitude,  m  is  the  zonal 

wave  number  out  to 

m  -  72  in 

this  case  (note  that  m  -  0  is  identically  zero 

because  the 

zonal  moan  was  removed).  Next,  the  normalized  zonal  power  '  (  ;  )  we. 


calculated  from  2  ,  2 

a  +  b 

< ; )  ,m  (3) 

m  2.2. 

El  a  *  o  ) 

,  m  m 
m=  1 

Thus  the  fractional  c  >nt  i  i but i ons  ovia  all  wave  numbers  lot  i  ;  at  te¬ 
rn 

ular  latitude  should  add  up  to  unity.  Figs.  4a  d  show  the  normalized 
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power  plotted  vs.  wave  number  for  both  the  Flattery  analysis  (,..'L3)  fields 

and  the  observation-corrected  fields  for  temperature  and  specific  humidity 

for  JUL.  This  spectrum  was  calculated  using  the  average  of  the  47.5°, 

50.0°,  and  52.5°  N  latitude  values,  from  which  the  departure  from  the  zonal 

average  was  calculated  and  a  and  b  were  obtained.  This  latitude  band 

m  m 

was  chosen  because  it  consistently  contained  the  greatest  number  of  ob¬ 
servations,  and  therefore  should  be  the  region  of  greatest  difference  be¬ 
tween  the  LVL3  field  and  the  corrected  field.  As  expected,  the  greatest 
difference  appears  at  wave  numbers  25  and  greater  especially  in  the  tem¬ 
perature,  because  of  the  truncation  of  the  Flattery  analysis  at  wave  num¬ 
ber  24.  The  successive  correction  procedure  puts  much  of  the  higher  wave 
number  variation  back  into  the  spectra  without  making  much  of  a  change 
at  lower  wave  numbers. 

Figs.  5a- f  illustrate  the  normalized  power  of  q  plotted  with  that  of 
T,  u,  and  v  for  comparison  at  this  same  latitude  band  for  JAN.  All  of 
the  plots  are  highly  oscillatory,  and  some  smoothing  of  this  phenomenon 
would  no  doubt  result  from  an  average  of  such  statistics  calculated  for 
a  group  of  observations.  The  linear  correlation  coefficient  value  r  is 
defined  by 


r 


X1X2  -  X1X2 


[  2  -2  2  -  2  "1  5 
[X!  -  X1  )(X2  '  X2  > 


(4) 


where  X  and  X.,  are  the  parameter  values  being  compared  (Fanofsky  and 
Frier,  lc)b8'*).  It  was  computed  using  the  normalized  power  for  the  two 
parameters  plotted  using  only  wave  numbers  one  through  12,  since  values 
for  greater  wave  numbers  approach  the  noise  level.  Such  a  value  gives 
some  basis  for  quantitative  comparison  between  T,  u,  and  v  as  to  which 
might  represent  Lite  best  approximation  to  the  spectral  composition  of  g. 

A  comparison  of  the  r  values  in  Table  2  shows  that  T  has  a  much  better 
correlation  with  g  than  do  u  and  v  in  all  cases.  Apparently,  another 
quantity  that  is  advected,  in  this  case  temperature,  is  more  closely  cor¬ 
related  to  specific  humidity  than  is  tlv  agent  responsible  for  the  ndveo- 
t i on ,  namely  u  and  v.  F'erhaps  the  greatest  contributing  factor  to  this 
relationship  is  the  tendency  ot  both  T  and  o  t<  h.nm  t ho  greatest  power 


9.  Fanofsky,  H.  A.,  and  ('..  W.  Brier,  19B8:  Some  Applications  of  Sta¬ 
tistics  to  Meteorology.  The  Pennsylvania  State  University,  Uni¬ 
versity  Park,  PA.  224  pp. 
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Fig.  4a.  Normalized  power  for  departures  from  zonal  mean  plotted 
as  a  function  of  wave  number  for  an  average  of  three  latitude  values 
centered  at  50  N.  Both  LVL3  and  corrected  field  results  are  plotted. 
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Fig.  4b.  Normalized  power  for  departures  from  zonal  mean  plotted 


as  a  function  of  wave  number  for  an  average  of  three  latitude  values 
centered  at  50°  N.  Both  LVL3  and  corrected  field  results  are  plotted 
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Fig.  4c.  Normalized  power  for  departures  from  zonal  mean  plotted 
as  a  function  of  wave  number  for  an  average  of  three  latitude  values 
centered  at  50  N.  Both  LVL3  and  corrected  field  results  are  plotted. 
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Fig.  4d.  Normalized  power  for  departures  from  zonal  moan  plotted 
as  a  function  of  wave  number  for  an  average  of  three  latitude  values 
centered  at  50°  N.  Both  LVL3  and  corrected  field  results  are  plotted. 
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Fig.  5a.  Normalized  power  for  departures  from  zonal  mean  for  a 
three  latitude  average  centered  at  50  N.  Values  for  corrected 
fields  only  are  plotted. 
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Fig.  5b.  Normalized  power  for  departures  from  zonal  mean  for  a 
three  latitude  average  centered  at  50°  N.  Values  for  corrected 
fields  only  are  plotted. 
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Fig.  5c.  Normalized  power  for  departures  from  zonal  mean  for  a 
three  latitude  average  centered  at  50  N.  Values  for  corrected 
fields  only  are  plotted. 
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Fig.  5d.  Normalized  power  for  departures  from  zonal  mean  for  a 
three  latitude  average  centered  at  50  N.  Values  for  corrected 
fields  only  are  plotted. 
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at  wave  number  2  and  then  decline  more  or  less  steadily  with  increasing 
wave  number.  The  greatest  amplitudes  in  the  u  and  v  spectra  tended  to 
be  at  somewhat  larger  wave  numbers  or  were  distributed  over  several  wave 
numbers.  This  seems  to  be  especially  true  for  v. 

TABLE  2.  CORRELATION  COEFFICIENT  VALUES  FOR  ZONAL  SPECTRA. 


e  Level 

Date 

q,  T 

q,  U 

q,  v 

850 

1/15/78 

.93 

.64 

.31 

7/20/78 

.97 

.42 

.12 

500 

1/15/78 

.86 

.59 

.49 

7/20/78 

.78 

.60 

.57 

The  normalized  zonal  power  spectrum  was  calculated  for  all  latitudes 
and  plotted  in  log^  form  on  latitude  vs.  wave  number  diagrams  as  seen  in 
Figs.  Ca-i.  The  contours  represent  the  fractional  contribution  of  each 
wave  number  to  the  total  power  at  each  latitude.  For  example,  a  point 
lying  to  the  left  of  the  -2  contour  indicates  that  that  wave  number  con¬ 
tributes  at  least  1%  of  the  total  power  at  a  particular  latitude.  In 
these  diagrams,  the  wave  numbers  correspond  to  the  number  of  waves  around 
each  latitude  circle  irrespective  of  convergence  of  the  meridians  near 
the  poles.  Therefore,  a  particular  wave  number  at  high  latitudes  will 
correspond  to  a  wave  of  shorter  length  than  the  same  wave  number  at  the 
equator  by  a  factor  of  cos  cf>.  This  is  reflected  in  the  fact  that  in  all 
diagrams  the  rate  at  which  the  power  falls  off  with  wave  number  increases 
toward  the  poles.  The  point  of  comparison  should  be  between  diagrams  at 
particular  choices  of  latitude,  and  not  between  latitudes  (except  for 
large  departures  from  the  cos  <ji  variation)  on  any  one  diaqram. 

In  comparing  the  four  q  diagrams,  we  see  no  appreciable  variation 
of  the  zonal  spectra  with  height;  similarly,  no  significant  variation  with 
season  was  observed  from  the  JUT.  plots  (not  shown).  Aside  from  the  differ¬ 
ences  in  latitudinal  location  of  the  spikes  in  the  -4  contour,  the  q  and  T 
diagrams  even  look  quite  similar.  The  spikes  indicate  latitude  bands 
where  the  variation  of  the  quantity  is  spread  out  over  a  greater  range  of 
wave  numbers  and  not  concentrated  ro  heavily  in  the  Iona  waves.  The  spike 
on  the  T  diagram  has  its  center  on  '  at  i*- tides  between  0°  and  20°  N,  whi  le 
the  q  diaqram  shows  two  lesser  spikes,  one  in  each  hemisphere  and  both 
of  which  seem  to  change  latitude  slightly  with  t  imo  and  height.  The  dia- 


Fig.  6b.  Zonal  power  spectra  in  powers  of  10  for  departures  from 
zonal  mean  for  500  mb  specific  humidity.  The  fractional  contribu 
tion  summed  over  m  for  each  latitude  should  equal  unity. 
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Fig.  6d .  Zonal  power  spectra  in  powers  of  10  for  departures  from 
zonal  mean  for  500  mb  temperature.  The  fractional  contribution 
summed  over  m  for  each  latitude  should  equal  unity. 
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Fig.  6f.  Zonal  power  spectra  in  powers  of  10  for  departures  from 
zonal  mean  for  850  mb  zonal  velocity.  The  fractional  contribution 
summed  over  m  for  each  latitude  should  equal  unity. 
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Fig.  6g.  Zonal  power  spectra  in  powers  of  10  for  departures  from 
zonal  mean  for  500  mb  zonal  velocity.  The  fractional  contribution 
summed  over  m  for  each  latitude  should  equal  unity. 
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Fiq.  6i.  Zonal  power  spectra  in  powers  of  10  for  departures  from 
zonal  mean  for  500  mb  meridional  velocity.  The  fractional  contri¬ 
bution  summed  over  m  for  each  latitude  should  equal  unity. 
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grams  for  Z  (500  mb  only)  show  a  greater  concentration  of  power  in  low 
wave  numbers  than  do  q  and  T.  A  comparison  of  the  u  diagrams  with  q 
shows  a  tendency  of  the  -2  contour  on  the  u  diagram  to  lie  somewhat  to 
the  left  of  that  of  the  q  diagram,  especially  at  middle  latitudes  in 
both  hemispheres.  This  is  clearly  not  the  case  with  the  v  spectra, 
which  in  some  cases  show  even  more  wave  numbers  included  in  the  1%  con¬ 
tribution  category  than  for  q.  From  these  qualitative  comparisons  it  can 
be  tentatively  concluded  that  temperature  seems  to  have  a  zonal  spectrum 
most  like  that  of  q,  where  u  tends  to  show  slightly  more  concentration 
of  power  in  low  wave  numbers  and  v  a  tendency  towards  greater  contribu¬ 
tions  from  higher  zonal  wave  numbers. 

The  total  power  of  the  zonal  eddies,  which  is  the  denominator  in  Eu. 
(3) ,  is  plotted  along  with  zonal  mean  power  on  a  log^  axis  as  a  function 
of  latitude  in  Figs.  7a- j.  The  fields  used  to  obtain  these  power  plots 
were  the  departures  from  the  global  mean  as  defined  in  Eqs.  (6)  and  (7). 
Total  eddy  power  for  q*  is  clearly  unique  in  its  latitudinal  variation 
when  compared  to  that  of  the  other  parameters.  The  power  remains  high 
in  a  broad  band  of  latitudes  and  falls  off  at  both  poles.  The  width  of 
this  band  of  high  power  is  greater  at  850  mb  than  at  500  mb  for  both 
dates.  At  both  levels,  the  total  eddy  power  is  greater  in  the  summer 
hemisphere  than  it  is  in  the  winter  hemisphere.  The  fact  that  this  pow¬ 
er  difference  between  the  seasons  is  greater  in  the  Northern  Hemisphere 
for  both  levels  is  probably  due  to  the  greater  number  of  observations  in 
the  Northern  Hemisphere. 

The  total  eddy  power  —  latitude  curves  for  T*  and  Z*  both  show  a 
bimodal  structure  witli  relative  minima  at  the  equator  and  poles  and  rela¬ 
tive  maxima  at  mid-latitudes.  Power  due  to  zonal  eddies  is  greater  in 
the  winter  Hemisphere,  as  one  would  expect  from  looking  at  weather  maps, 
but  as  such  it  is  just  opposite  from  the  moisture  variance.  As  in  the 
case  of  the  total  eddy  power  for  moisture,  the  winter  hemisphere-summer 
hemisphere  difference  is  greater  in  the  Northern  Hemisphere.  As  expected, 
the  total  power  curves  for  Z*  500  show  the  same  characteristics  as  the  T* 
curves . 

The  fact  that  the  total  power  of  the  zonal  eddy  specific  humidity 
appears  to  be  greater  in  the  summer  hemisphere  in  contrast  to  the  other 
quantities  may  possibly  be  explained  by  the  limiting  factor  of  saturation. 
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Q  850  1/15/78 

-  TOTAL  EDDY  POWER 

x-x-x  ZONAL  MEAN  POWER 


»*o 


Fig.  7a.  Total  power  due  to  the  departures  from  the  zon^l  mean  (de¬ 
nominator  in  Eq.  (3))  and  power  due  to  the  zonal  mean  (a,-,  (<p))  for  850 
mb  specific  humidity  for  JAN.  Quantities  are  calculated  from  the  de¬ 
partures  from  the  global  mean. 


I 


TOTAL  POWER  (G/KG)«*2 


Q  850  7/20/78 

- — —  -  TOTAL  EDDY  POWER 
x-x-x  ZONAL  MEAN  POWER 


«*o 


i.  7L.  Tot.i!  power  tue  to  the  departures  from  th^  zonal  mean  (denomin 
t  in  Ei | .  (  1))  and  power  due  to  the  zonal  mean  (a"  ($))  for  850  mb  spe 
:  i('  humidity  tor  JUL.  (quantities  are  calculated  from  the  departures 
»m  the  1  'b-i  1  mean. 
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Kig.  7c.  Total  powet  iue  to  the  departures  from  thq  zonal  mean  ( demon in- 
ator  in  Eq.  (3))  and  power  due  to  the  zonal  mean  (a"  (<(>))  for  500  mb  spe- 
cific  humidity  for  JAN.  Quantities  are  calculated  from  the  departures 
1  rum  the  g  1 1  jb.  1 1  mean  . 
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Fig.  7d.  Total  power  due  to  the  departures  from  th^  .^onal  mean  (denomin 
ator  in  Eq.  (3))  and  power  due  to  the  zonal  mean  (a  (<{>/)  for  500  mb  spe 
cific  humidity  for  JUL.  Quantities  are  calculated  from  the  departures 
from  the  global  mean. 
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Fig.  7e.  Total  power  due  to  the  departures  from  th^  zonal  mean  (denomin¬ 
ator  in  Eq.  (3))  and  power  due  to  the  zonal  mean  (a^  ( 4> ) )  for  850  mb 
temperature  for  JAN.  Quantities  are  calculated  from  the  departures  from 
the  global  mean. 
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Fig.  7f.  Total  power  due  to  the  departures  from  th^  zonal  mean  (denomin¬ 
ator  in  Eq.  (3))  and  power  due  to  the  zonal  mean  (a  ($))  for  850  mb 
temperature  for  JUL.  Quantities  are  calculated  from  the  departures  from 
the  global  mean. 
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Pig.  7g.  Total  power  due  to  the  departures  from  t^e  zonal  mean  (denomina¬ 
tor  in  Eq.  (3))  and  power  due  to  the  zonal  mean  (a  ( tf> ) )  for  500  mb  tem¬ 
perature  for  JAN.  Quantities  are  calculated  from  the  departures  from  the 
global  mean. 
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Fig.  7h.  Total  power  due  to  the  departures  from  t^e  zonal  mean  (denomina¬ 
tor  in  Eq.  (3))  and  power  due  to  the  zonal  mean  (a  (<j>))  for  500  mb  tem¬ 
perature  for  JUL.  Quantities  are  calculated  from  the  departures  from  the 
global  mean. 
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Fig.  7i.  Total  power  due  to  the  departures  from  t^e  zonal  mean  (denomina 
tor  in  Eq.  (3))  and  power  due  to  the  zonal  mean  (a  (<)>))  for  500  mb  geo¬ 
potential  height  in  JAN.  Quantities  are  calculated  from  the  departures 
from  the  global  mean. 
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Northern  Hemisphere  average  values,  given  in  Table  3,  show  that  while  tem¬ 
perature  increases  by  only  about  4  percent  between  winter  and  summer,  the 
saturation  specific  humidity  qg  (which  acts  as  the  upper  limit  for  q)  in¬ 
creases  by  73  percent  at  850  mb  and  59  percent  at  500  mb.  Thus,  the  max¬ 
imum  allowable  departure  from  the  average  is  about  2.6  g/kg  in  the  winter 
compared  to  4.9  g/kg  in  the  summer  at  850  mb.  The  hemispheric  average 
relative  humidity,  however,  for  the  two  cases  is  about  the  same,  suggest¬ 
ing  that  this  quantity  might  act  as  a  less  severe  limit  to  the  departure 
from  the  mean  in  the  winter  time.  To  investigate  this,  the  total  eddy 
power  curves  for  RH  based  on  the  ratio  of  q  to  qg  were  generated  and  ap¬ 
pear  in  Figs.  8a-d.  The  fact  that  the  winter-summer  hemisphere  differ¬ 
ence  is  lessened  for  RH  may  indicate  that  lower  temperatures  tend  to  re¬ 
duce  the  total  variance  about  the  zonal  mean  more  for  specific  humidity 
than  for  relative  humidity,  which  would  appear  to  explain  the  lower  var¬ 
iance  for  q  in  the  winter  hemisphere. 

TABLE  3.  NORTHERN  HEMISPHERE  AVERAGES. 


Level 

Date 

q(g  kg  1) 

T(°K) 

qg(g  kg  1) 

RH 

850 

1/15/78 

5.35952 

276.49600 

8.00103 

.66985 

850 

7/20/78 

8.95516 

288.18456 

13.83457 

.64730 

500 

1/15/78 

1.23598 

254.71750 

2.67001 

.46291 

500 

7/20/78 

2.13636 

267.28384 

4.25537 

.50204 

Following  the  zonal  harmonic  analysis,  a  surface  spherical  harmonic 
analysis  was  performed  on  the  fields.  This  was  carried  out  as  an  expan¬ 
sion  ot  the  departure  from  the  qlobal  mean  given  by 

*  N  M  m  m 

A  ( 0 ,  X )  =  £  XI  m  cos  mX  +  B  sin  mX)  (sin  <j>)  (5) 

n=|m|  m=0 


whe  re 


A  «t>,X)  =  A  «M)  -  [A] 

/TT/2  _ 

A  ($)  cos  <j>  d<J> 

h 


W  -  m  +  N 


total  power 


Fig.  8a.  Total  zonal  power  for  the  departures  from  the  zonal  average 
for  the  ratio  q/q^.  Both  q  and  T  used  in  the  calculations  are  the  cor¬ 
rected  field  values. 
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Fig.  8b.  Total  zonal  power  for  the  departures  from  the  zonal  average  for 
the  ratio  q/q  .  Both  q  and  T  used  in  the  calculations  are  the  corrected 
field  values. 
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Fig.  8c.  Total  zonal  power  for  the  departures  from  the  zonal  average 
for  the  ratio  q/qs.  Both  q  and  T  used  in  the  calculations  are  the  cor 
rected  field  values. 
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Fig.  8d.  Total  zonal  power  for  the  departures  from  the  zonal  average 
for  the  ratio  q/q  .  Both  q  and  T  used  in  the  calculations  are  the  cor¬ 
rected  field  values. 
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whore  N  represents  the  number  of  n  modes  and  is  a  constant  for  a  particu¬ 
lar  trunction,  i’^m  is  the  associated  Leqendre  function,  and  A  is  the  zonal 
mean.  The  expansion  in  Eq.  (5)  becomes  exact  for  M,  M  -*  but  the  number 
of  arid  | oints  around  latitude  and  lonqitude  circles  sets  a  maximum  limit 
of  72.  From  the  preceding  zonal  harmonic  analyses,  it  is  clear  that  amp¬ 
litude  approaches  the  noise  level  at  m  =  32;  so  a  truncation  limit  of 
M  =  32  was  chosen.  To  evaluate  the  integral  in  Eq.  (7) ,  the  global  mean 
of  the  qua  tity,  a  Gaussian  quadrature  integration  in  the  form 


Ak  Wk 


was  used,  where  are  the  corresponding  Gaussian  weights.  The  2.5°  lat¬ 
itude  values  were  interpolated  linearly  to  the  Gaussian  latitudes,  the 

■k 

quadrature  was  performed,  then  the  values  A  were  obtained  at  the  Gaussiar 

* 

latitudes  using  Eq.  (6).  Once  the  departure  values  A  were  obtained,  a 

Fast  Fourier  transform  was  used  to  evaluate  the  coefficients  Y  and 

mk 

5  of  the  expansion 
"k 

*  M 

A  (<t>,  ,A)  =  £  y  cos  mX  +  6  sin  mX  (10) 

m=0  mk  \ 

where  k  represents  the  indices  of  the  Gaussian  latitudes,  and  y  ,  6 

m  m 

are  given  by  x  K 

N 

Em  m  , 

a  P  (y,  ) 
m,  i  I  n  n  k 


6  P  (u.  ) 
n  n  k 


where  y,  =  sin  sk  -  Using  the  orthonormality  of  P  ,  we  have 


(U.  )  du 
k 


(u.  )  du. 

k 


08 


These  integrals  were  also  evaluated  using  the  Gaussian  quadrature  tech¬ 
nique.  For  M  =  32  and  64  Gaussian  latitudes,  the  highest  order  P  m  that 

32  n 

can  be  integrated  exactly  is  P,,  .  Therefore  M  =  63  and  N  =  31  were 

b  3 

chosen  to  complete  the  rhomboidal  truncation  of  the  expansion  in  Eq.  (5) . 
The  forms  of  the  quadrature  for  Eqs.  (12)  are 
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m 

a  =  E 

n 


,  .  Y  P  (y,  )  W 

k=l  n  k  Jt 
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6  m  =  E  S  Pm  (u  )  w.  . 
n  k=l  mk  n  k  k 


(13) 


The  normalized  power 


A  W 

A  was 
n 


then  evaluated  in  the  form 


,  m.  2  ,  .  nr  2 

(a  )  +  (8  ) 

n  n 


| m | +N  M 

E  E 

n=|m|  m=0 


r,  m.  2 

K  >  ♦ 


(8  V] 

n  J 


(14) 


Figs.  9a- i  display  the  results  of  the  spherical  harmonic  analyses 

for  JAN.  The  quantity  n-m  represents  the  number  of  times  a  particular 

harmonic  passes  through  zero  between  the  two  poles;  it  can  thus  be 

thought  of  as  a  sort  of  meridional  wave  number.  lor  this  reason  the 

logarithm  of  A  m  was  (lotted  on  n-m  vs.  m  diagrams.  As  in  the  zonal 
n 

spectrum  figures,  the  log  values  are  plotted,  so  that  the  -3  contour 
refers  to  0.1  percent  contribution  to  the  total  power.  As  in  the  case  of 
the  zonal  spectra,  there  is  very  little  variation  in  the  spherical  spec¬ 
tra  witli  season  or  height.  However,  when  comnarod  to  the  q  figures,  the 
temperature  spectra  show  a  greater  concentration  of  power  in  low  wave  num¬ 
bers  both  zonally  and  meridionally.  More  of  the  total  power  is  contributed 
by  higher  wave  numbers  in  q,  so  that  short  wave  length  variations  are  more 
important,  in  the  moisture  variation  than  in  the  temperature  variation. 

Tliis  difference  in  power  distribution  is  even  more  evident  in  the  Z  500 
figure,  especially  in  the  concentration  of  power  in  the  low  zonal  wave  num¬ 
bers  which  was  seen  in  the  zonal  spectra.  The  power  distribution  for  u 
is  characterized  by  a  marked  spreading  of  the  power  toward  higher  meridional 
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Fig.  9a.  Spherical  harmonic  analysis  for  the  departures  from  the 
global  average.  Contour  values  represent  powers  of  10  which  are 
the  fractional  contribution  to  the  total  power  by  each  wave  number 
combination . 
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Fig.  9c.  Spherical  harmonic  analysis  for  the  departures  from  the 
global  average.  Contour  values  represent  powers  of  10  which  are 
the  fractional  contribution  to  the  total  power  by  each  wave  number 
combination. 


Fig.  9d.  Spherical  harmonic  analysis  for  the  departures  from  the 
global  average.  Contour  values  represent  powers  of  10  which  are 
the  fractional  contribution  to  the  total  power  by  each  wave  numbe 
combination . 
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Fig.  Sie .  Spherical  harmonic  analysis  for  the  departures  from  the 
global  average.  Contour  values  represent  powers  of  ]0  which  are 
tile  fractional  contribution  to  the  total  power  by  each  wave  number 
combination . 
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Fig.  9g.  Spherical  harmonic  analysis  for  the  departures  from  the 
global  average.  Contour  values  represent  powers  of  10  which  are 
the  fractional  contribution  to  the  total  power  by  each  wave  number 
combination . 


wave  numbers  and  greater  concentration  of  power  in  low  zonal  wave  numbers. 
The  v  diagrams  indicate  a  somewhat  greater  spectral  spreading  in  both  dir¬ 
ections  than  q,  although  less  meridionally  than  u.  Also,  because  of  the 
very  small  magnitudes  of  the  zonal  average  of  v,  the  power  contribution 
from  the  m  =  0,  n  -  m  =  1,  2,  3,  4,  ...  wave  numbers  is  very  small. 

The  greater  concentration  of  power  in  long  waves  for  temperature  as 
compared  with  specific  humidity  is  shown  graphically  in  Figs.  jOa-d. 

Here,  the  cumulative  fraction  of  the  power  contributed  by  all  the  compon¬ 
ents  that  have  both  n-m  and  m  less  than  or  equal  to  the  value  of  "wave 
number"  on  the  abscissa  is  plotted  against  wave  number.  Thus,  for  example 
the  value  plotted  on  the  figure  at  wave  number  2  represents  the  sum  of  the 
fractional  power  contributions  i rom  all  wave  number  combinations  such  that 
both  m  and  n-m  are  less  than  or  equal  to  2.  The  difference  between  values 
for  consecutive  wave  numbers  indicates  the  fraction  contributed  by  the 
larger  wave  number.  Since  the  analysis  was  per  formed  on  departures  from 
the  global  average,  there  is  no  wave  number  zero  contribution. 

From  the  diagrams  it  is  evident  that  while  there  is  little  seasonal 
or  height  variation  in  the  power  distribution  for  temperature,  there  is 
a  distinct  variation  with  height  in  the  q  curves.  Wave  numbers  1  and  2 
together  contribute  some  60  -  65  percent  of  the  total  power  at  850  mb, 
but  less  than  40  percent  at  500  mb  for  both  dates.  In  all  cases,  it  is 
clear  that  a  greater  portion  of  power  is  accounted  for  in  the  first  two 
wave  numbers  for  temperature  than  for  q.  These  graphs  demonstrate  the 
error  incurred  ir  a  rhomboidal  truncation  at  a  given  wave  number.  For 
example,  choosing  a  wave  number  10  truncation  in  an  analysis  would  be 
quite  satisfactory  for  temperature  since  in  all  cases  well  over  90  per¬ 
cent  of  the  variance  would  be  accounted  for.  However,  this  would  not  be 
a  satisfactory  choice  for  q,  since  a  lesser  portion  of  the  variance  would 
be  included.  In  all  cases  for  both  q  and  T,  the  wave  number  2  contribu¬ 
tion  is  greatest.  The  cumulative  fractional  power  curves  for  u  and  v 
(Figs,  lla-d)  reveal  a  significantly  different  picture.  The  larger  contri¬ 
butions  are  at  higher  wave  numbers  and  no  single  wave  number  is  so  clearly 
dominant  as  is  the  case  in  the  q  and  T  curves. 

Because  the  fields  analyzed  spectrally  in  this  study  are  corrected 
Flattery  analyses,  it  is  conceivable  that  some  characteristics  of  the 
Flattery  analysis  procedure  may  be  evident  in  the  fields.  One  feature  of 
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Fig.  10a.  Cumulative  fractional  power  for  the  spherical  harmonic  analy¬ 
ses.  Wave  numbers  refer  to  equal  values  of  m  an<i  n-m,  and  fractional 
amount  plotted  refers  to  contribution  by  all  wave  numbers  less  than  or 
equal  to  the  indicated  wave  number 
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Fig.  10c.  Cumulative  fractional  power  for  the  spherical  harmonic 
analyses.  Wave  numbers  refer  to  equal  values  of  m  and  n-m,  and 
fractional  amount  plotted  refers  to  contribution  by  all  wave  num¬ 
bers  less  than  or  equal  to  the  indicated  wave  number. 
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Fig.  lOd.  Cumulative  fractional  power  for  the  spherical  harmonic 
analyses.  Wave  numbers  refer  to  equal  values  of  m  and  n-m,  and 
fractional  amount  plotted  refers  to  contribution  by  all  wave  num¬ 
bers  less  than  or  equal  to  the  indicated  wave  number. 
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Fig.  11a.  cumulative  fractional  power  for  the  spherical  harnnn 
analyses  of  the  wind  components. 
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Fig.  11c.  Cumulative  fractional  power  for  the  spherical  harmonic 
analyses  of  the  wind  components. 
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Fig.  lid.  Cumulative  fractional  power  for  the  spherical  harmonic 
analyses  of  the  wind  components. 
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Table  4  qivos  the  fractional  contributions  from  several  wave  number 


band. 5  for  the  four  quantities  q,  T,  u,  v.  For  the  lowest  wave  number 
band,  the  power  contribution  for  u  and  v  is  very  small  compared  to  T 
and  q  850,  an1  somewhat  smaller  than  the  q  500  cases.  However,  the  tables 
are  turned  in  the  3-10  wave  number  band,  where  the  contributions  for  u 
and  v  are  substantially  larqer.  The  contributions  for  q,  u,  and  v  are 
closer  in  the  highest  wave  number  band.  In  summary,  q  is  unique  in  show¬ 
ing  much  more  variation  with  height  than  the  other  three  parameters. 

There  are  no  great  changes  in  any  of  the  quantities  with  season,  and  q 
stands  quite  alone  in  its  distribution  of  power  with  wave  number,  lying 


somewhere  between  the  high  concentration  in  low 

numbers  for  T 

and  the 

greater  distribution 

of  power  from 

shorter  waves 

for 

the  wind 

components . 

TABLE  4.  FRACTIONAL 

CONTRIBUTIONS 

TO  TOTAL  POWER  IN 

SPHERICAL  HARMONICS 

SPECTRUM  FROM  SEVERAL  WAVE  NUMBER  BANDS. 

DATE  LEVEL 

WAVE  NUMBER 

q 

T 

u 

V 

BAND 

1-2 

.63 

.80 

.12 

.01 

850 

3-10 

.28 

.  16 

.70 

.67 

11-20 

.07 

.03 

.16 

.28 

1/15/78 

1-2 

.  38 

.83 

.06 

.02 

500 

3-10 

.43 

.  15 

.77 

.76 

11-20 

.16 

.02 

.15 

.19 

1-2 

.65 

.83 

.11 

.05 

850 

3-10 

.22 

.  13 

.66 

.70 

11-20 

.  10 

.04 

.20 

.21 

7/20/78 

1-2 

.37 

.87 

.  12 

.01 

500 

3-10 

.  38 

.  11 

.72 

.77 

11-20 

.21 

.02 

.  14 

.  19 

The  preceding  discussion 

lumped  zonal 

and 

meridional  wave  num- 

bers  together  in  discussing  spectral  composition 

.  It 

is  also 

of  interest 

to  determine  the  largest  contribution  from  a  single  wave  number  combina- 

tion.  Table  5  shows 

the  fractional 

contribution 

f  rom 

the  n  = 

2,  m  =  0 

r 


I 


wave  for  each  of  the  quantities.  In  addition,  the  maximum  values  for 
any  single  wave  number  combination  for  u  and  v  are  given.  For  q,  T,  and 
Z,  these  contributions  stood  out  as  by  far  the  largest  of  any  wave  num¬ 
ber  combination.  This  is  expected  since  it  represents  the  maximum  at 
the  equator  and  minimum  at  the  poles  displayed  by  the  zonal  mean  of  these 
quantities.  The  zonal  wind  component  shows  its  maximum  power  in  the  n  =  4, 
m  =  0  harmonic,  reflecting  the  dominance  of  the  zonal  mean  in  its  varia¬ 
tion  between  wester lies  at  mid- latitudes  and  easterlies  in  the  tropics. 
Because  ;  zonal  mean  is  so  small,  v  shows  its  largest  power  at  zonal 
.vo  number  4  for  JAN  and  6  for  JUL,  and  for  the  variation  passing  through 
two  zeros  between  the  poles. 


TABLE  5.  FRACTIONAL  CONTRIBUTION  OF  n 


2,  m  =  0  HARMONIC  TO  TOTAL  POWER. 


Level 

Date 

q 

T 

Z 

u 

V 

850 

1/15/78 

. 54238 

.66164 

— 

. 10771 ( . 11622) * 

. 00014 (. 06318) 

500 

1/15/78 

.  29926 

. 72589 

. 75830 

.04684 ( . 14380) 

.00036(. 06103) 

850 

7/20/78 

.40235 

. 52509 

— 

,05337(. 06181) 

.00084(. 04251) 

5 1)0 

7/20/78 

. 24650 

. 59242 

. 52723 

.05083 ( . 12871) 

.00001 ( . 10358) 

*  Values  in  parentheses  for  u,  v  represent  the  largest  contribution  from 
a  single  harmonic:  for  u,  in  all  cases  the  n  =  4,  m  =  0  harmonic;  for  v, 
the  n  -  m  =  2,  m=4  harmonic  for  1/15  and  the  n-m=2,m=6  harmonic 
for  7/20. 


The  correlation  coefficient  given  in  Eq.  (4)  was  evaluated  between  the 
fractional  spherical  harmonic  power  spectrum  values  of  g  and  those  of  T, 
u,  and  v.  The  results  are  given  in  Table  6.  The  large  values  for  T  rela¬ 
tive  to  those  of  .j  and  v  come  about  directly  from  the  fact  that  both  q 
and  T  have  relit  rely  large  fractional  contributions  at  common  harmonics 
(especially  u  ,  m  -  O)  compared  t  o  those  of  u  and  v.  However,  the 
Linger  of  drawing  conclusion-  on  the  basis  of  correlations  alon-  is 
seen  when  one  recalls  the  differences  between  power  distribution  of  q 
and  T  in  the  cumulative  fractional  power  curves.  What  is  surprising 
is  that  the  greater  difference  in  power  distribution  between  q  and  T 
at  500  mb  than  at  P.50  mb  is  t  indicated  at  all  in  the  correlation 
coefficient  values.  It  appears  that  the  two-dimensional  correlation 
coefficient  values  are  heavily  influenced  by  the  coincidence  of  large 
values  at  i  list  a  very  ;  ew  components. 
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TABLE  U. 


CORRELATION  OF  T,  u,  iUMU  v  WITH  q  FOR  SPHERICAL  HARMONIC  SPECTRA. 


Level 

350 
500 
850 
500 

For  completeness,  the  values  r‘or  total  :  c*v!  ‘  a  i.Li  :  L 
monic  analyses  are  given  in  Table  7.  Here,  ;  omf;  t<;  i 

lack  of  change  of  total  power  with  season  which  ;  ov:  : 

four  parameters.  In  all  four  cases,  the  greater  v.irl.u.  , 
Northern  hemisphere  winter.  This  may  result  fron  ‘.lie  •  •  : 


Dat  e 

T 

u 

1/15/78 

.99 

.62 

,  xj  / 

1/ 1 5/ 78 

.  97 

.  30 

.  0i 

7/20/73 

.  97 

.  4i 

.  05 

7/20/78 

•  92 

•  38 

.  Od 

observations  and  therefore  greater  corrections  tc  a  re!."  -  ■  .  )  >-m>-  >*■>: : 

Flattery-analyzed  field.  While  there  are  more  c :  a  i .  ■  :  it. 

Morthert.  Hemisphere  also,  this  may  be  part  iul  ly  1  two  • 

Lance  of  q  in  the  winter  hemisphere  as  seen  *> om  •  he  sen..:  sis  "  pj;: 
a lyses . 


TABLE 

7 .  total. 

POWER  for 

SPHERICAL 

harmonic;-  . 

Level 

Date 

q  (C,/KO)  “ 

o  2 

T  (  K) 

Z  fM)  " 

;  M  »  " 

v  !M  " 

O-Q 

1/15/73 

41 .26366 

358. 15657 

-•  - 

' !  *  r,  f • 

12--.  -.i-i 

500 

1/15/78 

6. 17996 

297 . 63 16] 

1 8 1255.2 84  9 i 

5 3  ‘  . 9 "34  3 

291.941  '8 

850 

7/20/78 

43.22954 

249.25719 

— 

1-1 V .  i  ’>5  04 

J07.41.li-  • 

500 

7/20/78 

3.77466 

333.86061 

1  27  382.0069  7 

78  9.  12  7 

2 1  ”  .  4 ‘  ' 

IV.  Summary  and  Conclusions 

The  following  is  a  summary  of  the  similarities  and  differences  noted 
when  comparing  observed  and  analyzed  global  specific  humidity  fields  with 
those  of  temperature,  zonal  wind  component,  and  meridional  wind  component 

1.  In  physical  space  on  the  synoptic  scale  there  is  qood  correlation  be¬ 
tween  the  q  perturbations  and  T,  u,  and  v  perturbations.  Warm,  moist  air 
is  transported  poleward  and  cold,  dry  air  equatorward  by  the  larqe-scale 
eddies . 

2.  The  zonal  harmonic  spectra  for  q  are  more  similar  to  those  of  T  than 
they  are  for  u  and  v  in  that  most  of  the  power  is  contributed  by  low  wave 
number  (large  scale)  waves.  The  zonal  wind  appears  to  be  even  more  re- 
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stncted  to  lonq  waves  in  its  power  distribution,  while  the  v  power  seems 
to  be  more  spread  out  over  shorter  waves. 

> .  While  the  total  power  of  the  departures  from  the  zonal  mean  for  q 
exhibits  a  broaa  maximum  in  a  band  of  latitudes  about  the  equator,  the  '1 
total  power  shows  a  minimum  at  the  equator  and  the  poles  and  relative  max¬ 
ima  at  mid-iat i tides .  donal  eddy  power  curves  for  u  and  v  are  decidedly 
more  monotonic  with  latitude.  In  addition,  the  total  zonal  perturbation 
rowei  fot  is  somewhat  less  in  the  winter  hemisphere  than  in  trie  summer 
hemisphei  •.  ,  m  direct  contrast  to  T,  u,  and  v.  However,  total  eddy  power 
Tor  relative  humidity  in  the  two  hemispheres  is  about  the  sane. 

4.  The  spherical  harmonic  spectra  for  q  are  similar  to  those  for  temper¬ 
ature  in  havinq  most  of  their  power  contributed  by  the  zonal  mean  that 
has  a  sinqle  wave  lenqth  between  the  poles  (minima  at  the  poles,  maxima 
at  the  equator) .  The  u,  v  spectra  showed  maximum  power  contributed  1// 
snort  or  waves  in  both  directions. 

j.  However,  t  he  spherical  harmonic:  power  spectrum  ioi  ■  i  is  less  con con¬ 
i'  rated  in  low  wave  numbers  than  is  that  :  •  >r  T,  but  is  not  as  a:  lvad  out 
m  hi  ihe.r  wave  namoer  iai  either  direction  as  are  a  as  :  v. 
a.  for  a  sinalo  wave  number  represent inu  an  equal  numoei  or  zonal  waves 
..is  zeros  between  tiie  poles,  the  :  ectra  show  a  *  iani:  ieantlv  creator 

.•oiii.vntrarioi!  of  power  in  the  l-l  wave  number  bund  at.  .hit1  mb  •  hau  at  50- 
mb.  iiiis  neiuui  difference  was  not  eviaent  in  :  ,  a,  or  v.  hobo  c  ;•  the 
isantities  snows  much  variation  in  too  i  -owe  r  distribut  ion  cits,  season. 

in  qonoral ,  t.hu  tcim.orature  and  joes  etcnt.iaJ  no  l  i:it  vo!  i  .•  have 
t  ne  i  r  lower  too  hiahiy  concent  rated  in  the  iow  wavs  t.nmix  :  t.<  •  or  sc',; 


!'i  ' 

;.  veson  ta  t  u 

Ills  !<'*•  1 

! .  riio 

wind  component  s  seem  to 

have  too  much 

S  'OWt  i 

1  £  1 

in  til  wave 

numbers 

t  o  be  r 

•0| 1  resen  ta  t.  we  r  o  r  .; .  A 

.1  l\: .J  !  t  ,  it 

Si  1  ■«.. 

L  r. 

at  neither 

the  IkiS: 

i  field 

nor  tile  wind  component  s 

lh.il  VSCS  will 

a:  sa- 

ruLoiy  |  or  tray  tile  moisture  field. 

Tilt!  results  ol  tin  s  study  seep  ro  indicate  *  hot  t  :c.  .u.iii  i;« u  ;s  a  i « 
analysis  •ciicme  is,  roouirod  m  ordei  to  aenor  it  e  a  sat  i.s  .a  •  :  •  ,.r, 

i  e|  to.seiit.il,  ion  tor  i  toroeust  model  in)  at  .  fC'-isse  id  tin  i".i  ;  .<:■  .a  1  ■■  - 
t.rai  character  oi  sreci:  i:  humid  i  t.v,  i  ;  vik  mid'  i  .■  ;  i  .- 

variate  analysis  for  moisture  usinq  anv  .  »•  tin  tin  i  j  .u  aniotoi  s 
sidered  in  this  study. 


B.  THE  GENERATION  OF  HARMONICS  FOR  ANALYSIS  OF  METEOROLOGICAL  DATA  IN  THE 


VERTICAL  AND  HORIZONTAL 

To  prepare  meteoroloqical  observations  as  initial  conditions  for  the 
proposed  AFGL  spectral  model ,  sets  of  orthogonal  functions  were  qenerated 
to  act  as  bases  in  the  analysis  and  initialization  of  the  data.  The  func¬ 
tions  are  solutions  of  the  Laplace  tidal  equations  which  are  linearized 
versions  of  the  equations  governing  the  behavior  of  a  r  Luid  in  a  basic 
state  of  rest  but  subject  to  coriolis  forces  and  linear  pressure  gradients 
The  functions  used  in  the  meridional  and  zonal  directions  are  the  Hough 
harmonics  which  yield  Fourier  components  zonally  and  Hough  functions  meri¬ 
dionally.  The  vertical  direction  is  handled  by  a  set  of  empirical  ortho¬ 
gonal  functions  (EOF)  accord inn  to  the  method  oi  Holms trom  (1063). 

The  Hough  functions  are  derived  according  to  the  method  discussed  by 
Kasahara  ( ly/o,  1  V78)  *.  Three  functions  are  derinc-d  re  latino  to  the  velo¬ 
city  potential,  stie.im  function,  and  •leopotent  iaJ  hoi  hit:,  respectively: 

V'  - 

a.  ;  -  l  A  i;  (.. ) 

i-.  n  s  n  n 

i . .  y  =  if  ’  (  1  ) 

c .  /.  _  c  i  (  ) 

s  n  s  n  n 

wnore  i  =  »-i  and  I  are  the  Legendre  functions  o r  ..  sin  ,  wnoro  v  is 
n 

latitude,  s  is  ei’.o  zonal  index  and  n  is  the  mor  iai  :,a  2  in  ivy  using 

(i)  in  the  basic  o  piations  o:  mot  i  ->n ,  one  can  derive  r<  ssions  :  or  the 

set  oi  coott  icient  s{.<\  J,  «*roi  J..  these  solutions  divide  into  two 

categories.  —  symmetric  and  an  l  l -svmmet  r  i  c  .  '1  he  svmmetric  mod.-  .  :.*  ..  Is 

on  i  v  even  values  oi  u  -  s  tor  A  and  and  add  valuer,  t  or  B  .  'hi  ai.ti- 

1 1  : :  a 


11.  K  i  rah,  oi,  a.,  Normal  nnilf’K  of  nl  t  r<t  long  waves  in  tin:  at  most  here. 

Mon.  We  a.  KOv.  ,  1 U4  ,  i.py-ngn. 

,  l!j7d:  F irthei  -dud:  i*  on  a  spool  ril  modi  1  o*  the  liobai 
bai  otropi  o  j  riniit  ive  equat  ions  with  Hough  luiimonic  expansions,  a.  At  mo: 
s,  i.,  in,  .  is  j  . 

'M 


!  I-  '■!  ?■.  1  i  ,  ■  i  i  '  •  i  II'.: 


when  ‘''n  ■:  i  if:.’.  ■ .  •  t 1  '  so  J  : : :  ions  soinii'f  .  Her.-o  .  a  Ui  ffr  r  :  1 1  at  - 


;  ■••jch  is  unde  rial-tr  r  whereby  t  at-  '•  i  -iral  o  :a.it  i  ::  ■-  ivc  i  •.■wri  tt.cn  ‘or  s  - 

f  and  tiie  non-di  menr-iuna  1  p>e  rod,  C  ••  u .  t’ndor  these  rord.it  ion*- 

wiii  In  i’  is  governed  by  a  qeosti  oplii  c  rel  at  i  onsh.i :  .  The  goo f r i eient s  Bc 

n 

and  i’-  arc  then  relited  by  substituting  into  the  induced  e  stations.  Tho 
n 

o  a  o 

ettSN  i  ng  rolat  ionshii-  involves  only  each  B.  wit  it  i\  .  .ma  T.  lojvtna 

K  K  *-  J  K  -  1 

the  other  coefficients  arbitrary.  Their  selection  is  made  so  as  to  von- 

o 

der  the  entire  set  orthonormal.  By  choosing  a  s|  >ei:  5  f  i  c  i,  and  ill 

others  zero,  a  set  of  N  vectors  can  lie  generated,  then  or t  iionorma  1  ire. i  to 

/V  j ?  A  f  "o  .  ‘O 

qive  U  '  and  7.  .  The  lowest  order  qives  U  =  0  tor  all  while  7,  must 

o  o  o  o 

become  1//2  for  all  p  in  order  to  satisfy  orthonormality.  Solution  for 
the  non-rota tional  modes  (i.e.,  the  gravitational  modes)  is  straighttor- 
ward.  The  resulting  expressions  for  and  can  be  simplified  to  give 
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-t  S  £  B°P 1 
n=l  n  n 


V  -  ^  53  A0?'1 

n=l 

where  the  first  order  associated  Legendre  polynomial  is  introduced  because 
o  2  1 

of  the  identity  DP  (i-U  )  5  =  P  . 

n  n 

^  o  ^  ■> 

The  roles  present  special  problems  for  the  computation  of  L  -  and  V' 

.  s  s 


because  Eg.  (2)  seems  to  indicate  a  singularity  when 


At  close  in¬ 


spection,  however  ,  Eu.  (2)  presents  no  problems  when  s  •  i,  because  the 

2 

Legend re  : unctions  starting  with  P  are  polynomials  at  least  of  order  2 

2  -t 

in  u ,  while  the  denominator,  (1-U  )  *,  is  only  first  order  in  U.  Since 

pS(tl)  -  0  for  all  s  >  0,  there  is  no  Question  that  U?(il)  =  V^(±l)  =  0 
n  s  s 

tor  all  s  1.  When  s  =  1,  the  formula 

2  2  -h  1  o 

P  -  2u(l-U  )  P  +  n  (n  +  1)  P  ^  0 

n  n  n 

1  2  -h 

mav  be  used  to  derive  the  limit  of  P  (1-M  )  as  u-*-±  1  and  hence  the  values 

of  U  (il) ,  V  (11) ,  which  will  not  be  zero  at  the  poles.  For  s  =  0,  Eg. 

1  'I--'  ?  1 

(4)  indicates  that  U  (11)  =  V  (13)  =  P  (11)  -  0,  while  Eq.  (1)  indicates 
,  o  o  n 

that  ZK  ('])  i  0,  because  P°  ( t 1 )  4  0. 
o  n 

Truncation  of  the  sums  for  the  generation  of  the  Houqh  functions  was 
always  at  s  +  39.  The  number  of  modes  calculated  was  40:  25  rotational 
modes,  8  westward-travellinq  gravity  modes  and  7  eastward-travelling  modes. 
The  functions  can  be  calculated  for  any  latitude  and  stored  on  disc  or  tape 
for  use  in  the  analysis  of  data.  To  date,  the  Hough  functions  have  been 
generated  for  every  2.5°  latitude  and  for  64  Gaussian  latitudes.  The  lat¬ 
ter  have  been  used  to  test  the  orthoncrmal ity  of  the  Hough  vectors,  which 

-9 

has  been  verified  to  an  order  of  approximately  10 

The  2.5°  latitude  values  have  also  been  used  in  a  computer  program 
which  was  developed  to  take  wind  and  height  values  on  a  2.5  latitude-lon¬ 
gitude  grid  at  any  number  of  pressure  surfaces  in  the  vertical  and  calcu- 

s  s 

late  the  coefficients  a.  ,  and  in  the  expansions 

X/  f  K  X-  /  K 


24 

40 

12 

u  u,e,p) 

=  E 

E 

E 

s=0 

£=1 

k=  1 

24 

40 

12 

v  ( A  ,  0 ,  p) 

=  E 

E 

E 

s=0 

£=1 

k=l 

^  S  S 

E  [ae  k  cos  (sX)  4  k  sin  (sX)J  us  (sine>  fv(p) 


sin  (sA)  -  cos  (sA)]  vg  (sin0 )  F  (p) 


24  40  12 

Z  .  9  ,p)  =  £  2  £  [ 

s=0  H= 1  k=l 


S  S  T  V 

af  k  cos  (v-A)  +  i-:£  ^  sin  (s>.)J  (sin 0)  f'k(p) 


in  which  u,  v,  and  z  are  the  non-dimensionalized  wind  and  height  values 
related  to  the  actual  grid  point  values  u,  v,  h  by 

u  =  u/C  ,  v  =  v/C  ,  z  =  h/h 
g  g  o 

S 

where  =  (ghQ)  and  hQ  is  the  equivalent  height.  A  value  ol  10  km  was 

used  for  hQ.  The  F^(p)  are  the  EOF's  which  will  be  discussed  shortly.  be- 
s  s 

tails  on  how  ct^  ^  and  3^  k  were  obtained  are  given  in  Kasahara  (1970)  and 
his  method  was  followed  exactly  except  the  s  -  0  case  was  included,  as  was 
the  vertical  dimension.  Once  the  coefficient  values  were  obtained,  they 
were  used  in  the  expansions  Eq.  (5)  to  retrieve  the  values  u,  v,  and  z  and 
subsequently  u,  v  ,  h.  This  analysis  procedure  was  :  er formed  on  FGGE 
Level  III  (analyzed)  data  for  July  20,  1978  at  00  GMT.  The  results  showed 
qood  agreement  with  the  Level  III  data;  however,  no  examination  of  tile 
nature  of  the  difference  in  values  between  the  Level  III  analysis  and  our 
analysis  has  yet  been  conducted. 

Assimilations  with  the  Hough  functions  are  being  tested.  Some  tests 
have  already  been  made  using  known  analytic  functions  of  u.  Initial  re¬ 
sults  are  disappointing.  The  fact  that  Hough  functions  are  orthogonal 
does  not  guarantee  that  they  will  necessarily  defane  or  represent  an  arbi¬ 
trary  function.  Obviously,  functions  which  are  not  identically  zero  at 
the  poles  will  not  be  accurately  portrayed  by  Hough  functions  of  higher 
order.  In  fact,  there  are  probably  many  classes  of  functions  which  cannot 
be  represented  by  even  the  associated  Legendre  functions  of  higher  order. 
Even  for  those  functions  which  can  be  represented,  truncation  can  become  a 
serious  matter,  especially  for  higher  values  of  zonal  wavenumber  s. 

The  functions  used  in  the  tests  included  polynomials  in  p  of  various 
order,  exponential  functions  of  g,  and  sines  and  cosines  of  2ng.  Some  were 
forced  to  equal  zero  at  the  [>oles,  while  some  were  not.  Several  were  sym¬ 
metric  with  respect  to  the  equator,  while  others  were  antisymmetric  and  yet 
others  were  asymmetric.  When  these  functions  were  expanded  in  associated 
Legendre  functions,  tiic  representation  grew  worse  as  s  increased.  Also, 

results  were  inevitably  worse  for  odd  s  than  for  even  s,  probably  due  to 
s 

the  fact  that  P  is  a  pure  polynomial  in  p  if  s  is  even  but  contains  a  fac- 
n 
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tor  of  (l-p*-)  5  if'  s  is  odd.  When  expanding  in  Hough  harmonics,  similar 

test  functions  were  employed  for  u  (fi)  ,  v  (jj)  ,  and  h  (h)  and  the  resulting 

expansions  compared  with  the  original  values.  The  errors,  in  general, 

were  at  least  of  the  same  order  as  the  expansions  in  associated  Legendre 

polynomials  with  odd  s,  and  usually  even  higher.  With  variables  normalized 

-2  -1 

to  1,  errors  were  of  the  order  10  -  10  .  Thus  expansions  in  Hough  func¬ 

tions  could  engender  severe  errors  in  portraying  "true"  meteorological 
fields. 

As  an  example  of  Hough  analyses.  Figs.  12a-d  and  lja-d  show  • •. 
height  fields  and  their  approximation  by  Hough,  series  for  a  symnu.-.: 
and  anti -symmetric  function,  respectively.  The  symmetric  .’unction  was 
chosen  to  give  h  =  ip  -4p  *2,,  -5gN  i,  u  -e  (24li  -24u  +8b"'-10u)  r  C;i)  , 


and  v  -■  -4p  +p ,  where 


r  (P)  ~ 


2  >.  -  1 

(1-U  )  *  H  ,  !bj>0. 3 


'  ,'t .3,  U|-0.3 


The  representat  ion  for  u  was  taken  a:,  a  geostio;  hie  relationship  for  extra- 
tropical  regions.  Notice  that  the  functions  satisfy  the  requirement  that 
they  be  zero  at  the  poles.  Pigs.  12a,  b,  c,  and  d  are  for  s  =  0,  3,  5, 
and  6,  respectively.  Agreement  is  acceptable,  with  decreasing  accuracy 
for  increasing  s.  Regions  most  affected  are  neai  the  poles  and  the  equa¬ 
tor.  The  anti-symmetric  f unci,  tons  are  given  by : 

Z  _ 1 ,  2 

h  -  sin'  2/i pj  +  4  sin  2-rp,  u  =•  ‘  (6ir  sin*'  2rr(.  cos  2 up  +  8r  cos  2tp>  r  (u)  , 

v  =  *j  cos  2fTp  -  ‘j. 

Hough  expansions  are  shown  in  Figs.  13a,  b,  c,  ana  d  for  s  -  0,  3,  5,  and  t , 
respectively.  Here,  the  errors  are  larger  than  for  the  previous  case  and 
the  deterioration  with  increasing  s  is  quite  noticeable.  In  fact,  besides 
an  increase  in  error  there  is  also  a  phase  shift,  with  the  analyzed  ridges 
being  displaced  to  the  right  and  the  troughs  to  the  left.  The  reason  for 
these  errors  is  unclear.  It  could  be  a  factor  of  truncation,  inaccurate 
integration,  or  disansociat ion  of  the  windfield  from  the  height  field.  Rep¬ 
resentation  of  the  same  h  functions  by  Legendre  functions  loads  to  much  closer 
approximation  even  for  higher  s  values.  This  would  seem  to  indict  not  so 
much  the  truncation,  which  was  the  same  for  both  the  Legendre  and  Hough 
functions,  but  the  difficulty  in  balancing  waves  of  different  frequencies 
with  the  Hough  functions. 

The  vertical  harmonics,  based  on  KOFs,  are  calculated  separately  for 
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VRLUES 


of  the  function  for 


of  the  function 


Symmetric  function  versus  latitude  and  Hough  analysis 
:tion  for  S  =  5. 


FUNCTIONS  -  H- VALUES 


of  the  function  for 


An ti -symme tr ic:  function  versus  latitude 


each  set  of  data.  Given  any  set  of  data  defining  Z  (x,y,p) ,  where  x  and  y 
are  the  horizontal  variables  and  p  is  pressure,  it  is  possible  to  obtain 
a  set  of  EOFs,  F^ (p)  ,  using  the  calculus  of  variations.  The  first  order 
EOF,  F  (p) ,  is  simply  the  average  over  x  and  y  for  a  particular  p.  That  is, 


Fo(^  =  s 


Z(x,y,p)  dS, 


where  S  ’’-dicates  the  surface  area  over  the  entire  field  of  x  and  y.  It  is 


now  nee. 


ary  to  consider  only  the  departure  from  the  mean,  z(x,y,p)  = 


Z  ( x , y , p )  -  F'  (p)  .  By  expanding  z  such  that 


z(x,  y,  p)  =  £  (x»y)  F’  (p)  +  r  (x ,  y,  p)  , 

k=i  k  k 

and  minimizing  the  residual  r,  we  arrive  at  conditions  governing  and  F^. 
By  requiring  that  these  conditions  be  applied  sequentially,  and  that  F  be 

r  2 

normalized  so  that  /  F^  dp  =  1,  wc  get  the  following  system  ot  equations, 


z  F1  dp  = 


v  /-.  './ 


z2  dS 


J  (Z-Zi  Fj)  F2 


J  <z'zl  F!>  z2  dS  =  F2  JZ 2 


l/p  /  (2  -  E  \  Fk)  Fn  dp  =2,.  f  (2  -  £  VV  Zn  dS 


k  k  n 


F  f z2  dS. 
n  I  n 

•'S 


The  sequential  nature  of  the  calculation  guarantees  that  the  low  order  value: 

z  and  F^  carry  the  major  portion  of  the  variable  z(x,y,p)  and  higher- 

values  of  z  and  F  will  eventually  become  minimal.  The  Eqs.  (6)  also 
k  k 

guarantee  that  the  sets  J  and  {z^_  J  are  orthogonal  in  the  sense  that 


.  -Maul.*--  • 


1 

p 


F 


k 


F 

3 


dp 


<5?  and 

K 


z  ds 


»  0  if  k  =f  j. 


The  components  {z^}  and  {f^}  can  he  theoretically  solved  usinq  Eq.  (6), 

since  for  each  pair  of  unknowns  F  and  z  there  are  two  eouations.  Prac- 

k  k 

tically,  the  equations  can  only  be  solved  through  iteration  once  an  appro¬ 
priate  quadrature  is  chosen  to  represent  the  integrals.  For  the  surface 
into  :ral  a  simile  sum  of  all  points  is  substituted  as  an  approximation 
to  the  double  integral.  This  means  that  distribution  plays  no  rc lo  in  *  he 
vim:  at.it  ion.  Areas  which  are  data  void  will  receive  very  little  wei  us  i:. 
r  hi  *5  sort  of  approximation.  The  vertical  integration  is  approximated  using 
*  :.e  t:  gezoidal  rule  so  that 


/  ' 


2P 


J-l 

£ 

3=1 


(Fn(Ppl'  +  !rn(hn  ( P  j  + 1 


‘V* 


The  actual  procedure,  then,  is  to  firs,  calculate  Z  (x,  y ,  p)  at  all  man¬ 
datory  pressure  levels  by  taking  all  observations  of  Z  (x,  y,  p) ,  deter¬ 


mining-  Fo(pkl  =  ~ 


iiuess  for  F 


M 

6  2 

<■  '  is  then 


(xi'vY 


'k'  ' 


and  subtracting  F 


■‘rom  Z.  An  initial 


evaluated.  In  j  racticc,  the  expression 


„  ( 1 )  ,  , 
'i  (|r' 


1  - 


t 


was  used,  where  p  is  the  pressure  at  the  lowest  level  and  p  is  the  p res¬ 
ts  ‘  t 

sure  at  the  topmost  level.  Tins  first  guess  was  chosen  only  on  the  basis 

of  varying  monotonies  11 y  with  height;  no  consideration  was  given  to  its 

orthonormality  since  the  subsequent  iterations  for  F  are  always  normal - 

1  (1) 

ized.  The  first  integral  in  Eq.  (G)  is  now  approximates:  to  give  , 

where  the  superscript  indicates  thg  iteration  step.  The  inteoral 


/. 


d$  is  approximated  by 


ID  ,c 

/.  z  dS 


M 

L 

i  =  i 


,xryi'|,k) 


(2) 

yield  Fj  ' (p  ) 


for  each  k. 


and  similarly 
The  ratio  of  these  will 


dp  is  evaluated  by  the  trape¬ 


zoidal  rule  and  each  F  j  (p  )  divided  by  q2  in  order  to  achieve  the  rc- 

1  K  (2) 

quired  normalization.  The  new  F,  (p^)  are  now  re-inserted  into  the  first 
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integral  and  the  process  repeated.  Convergence  is  usually  very  swift  de¬ 
pending,  of  course,  on  the  convergence  criterion. 

The  drawback  with  this  method  is  that  the  vertical  modes  must  be  re¬ 
calculated  for  each  data  set.  Also,  there  is  little  theoretical  support 

for  the  calculation  of  EOFs,  but  alternate  methods  such  as  the  numerical 

1 2 

methods  recommended  by  Kasahara  (1981)  may  not  be  worth  the  effort. 

One  possible  alternative  may  be  the  eulcul  •‘ion  of  vertical  EOl's  which 

are  function  oi  the  mandatory  levels  only  and  need  not  be  updated  at 

ouch  ar.alys  .  time.  Such  a  method  was  attempted  by  Dixon  et  al.  (1972) 1 3 

1  -4 

tor  all  three  dimensions  with  limited  success 

The  status  of  the  current  analysis  method  may  thus  be  classified  as 
encouraging  but  subject  to  cautious  evaluation.  The  liouqh  analysis  does 
not  represent  an  arbitrary  field  of  variables  with  great  accuracy.  but 
Hough  functions  are  not  chosen  as  bases  because  of  their  accuracy.  Their 
primary  use  is  in  damping  the  growth  of  unwanted  gravity  waves.  The  ques¬ 
tion  is  whether  the  loss  of  accuracy  in  the  initial  conditions  will  depre¬ 
ciate  the  forecast  more  than  the  gravity  noise.  Also,  a  problem  arises 
when  interpolating  data  onto  a  fixed  grid  system  in  order  to  prepare  an 
analysis.  But.  tiiis  is  a  problem  related  to  data  distribution  and  is  corn- 


id.  Kasahara,  A.  and  K.  Puri ,  1981:  spectral  representation  of  three-di¬ 
mensional  global  data  by  expansion  in  normal  mode  functions.  Mon. 

Wo  a  .  Kov .  ,  1 i'9,  17-51 . 

13.  Dixon,  E.,  A.  Spackman ,  .1.  Jones,  and  A.  Francis,  1972:  1'ile  global 

analvsis  of  meteorological  data  using  orthogonal  j olynomial  base 
functions.  J_._  Atmos.  Sci  ■  ,  29,  609-022. 

14.  A  comparison  of  EOF  generation  t-.inq  fixed  levels  with  a  sequential 

derivation  of  coefficients  of  a  given  function,  similar  to  Holmstrom's 
method,  is  presented  in  a  report  by  Halberstam  ("An  Investigation  of 
Three  Methods  for  Spectral  Representation  of  Randomly  Distriouted 
Data,"  Scientific  Report  No.  1,  1981,  issued  under  this  contract). 
Results  show  that  although  EOF  generally  produce  more  accurate  approx¬ 
imations  than  sequential  derivation,  they  are  very  sensitive  to  spa¬ 
tial  distribution  of  data  points.  In  fact,  results  from  EOF  were  i- 
dentical  to  those  produced  from  a  straiqht forward  matrix  inversion 
where  the  known  quantities  were  given  at  the  observation  levels.  A 
matrix  of  this  form  is  known  to  be  ill-conditioned  because  of  its 
dependency  on  spatial  distribution,  and  when  the  matrix  solution 
fails,  the  EOF  method  will  fail  also. 
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mon  to  all  analysis  schemes.  The  vertical  analysis  has  been  used  opera¬ 
tionally  for  some  time  but  still  has  shortcominqs .  Despite  its  quick 
convenience,  the  scheme  requires  updatinq  at  each  analysis  time.  This 
could  become  a  costly  procedure.  It  also  fails  to  weiqh  data  from  data- 
poor  areas  more  heavily  than  data  from  data-ricn  areas.  some  of  these 
problems  can  be  overcome  by  adoptinq  certain  features  from  other  analysis 
schemes  such  as  optimal  interpolation.  Future  research  will  be  (eared 
to  investigating  techniques  used  by  other  analysis  procedures  to  overcome 
some  of  the  current  problems. 


i 
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SECTION  2. 

MESOSCALE  FORECASTING 


E 

14 
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A.  MESOSCALE  FORECAST  EXPERIMENT 


I.  Introduction 

The  need  for  developing  interactive  aids  for  formulating  short-range 
forecasts  led  to  SASC  participation  in  the  AFGL  Mesoscale  Forecast  Exper¬ 
iment  (MFE) .  Development  of  a  menu-driven  Mesoscale  Forecast  Facility 

1  2 
(MFF)  and  the  forecast  experiment  conducted  with  it  constituted  Phase  I 

of  the  MFE.  This  section  reports  on  the  follow-on  program  which  contin¬ 
ued  research  into  the  kinds  of  data,  displays,  and  capabilities  to  be 
supplied  to  interactive  facilities  (MFE  Phase  II) .  More  rigorous  fore¬ 
casting  sessions  lasting  an  entire  working  day  were  conducted  and  only 
two  predictands,  wind  and  cloud  amounts,  were  forecast.  Participation 
was  limited  to  six  forecasters.  Results  of  the  experiment  were  expected 
to  lead  to  development  of  techniques  and  software  to  enhance  interactive 
forecast  facilities. 

II.  Software  Modifications 

Software  changes  were  required  to  modify  the  MFF  (Phase  I)  program 
prior  to  the  start  of  Phase  II  forecasting  activities.  The  data  file 
used  to  store  forecast  and  skill  score  information,  FCSFIL,  had  to  be 
restructured  to  accommodate  changes  in  the  numbers  of  forecasters,  fore¬ 
cast  periods,  and  predictands.  Most  of  the  MFF  menu  was  eliminated  in 
compliance  with  the  Phase  I  conclusion  that  displays  of  conventional  and 
satellite  data  could  be  generated  more  easily  using  regular  McIDAS  com¬ 
mands.  The  chief  forecaster  position  was  eliminated,  leading  to  dele¬ 
tion  of  another  branch  of  the  menu.  Predictand  selection  and  review  was 
also  eliminated  in  the  Phase  II  format.  All  that  remained  of  the  MFF  menu 
were  the  forecast  entry  and  interrogation  modules,  both  of  which  required 
rewriting  because  of  changes  in  the  file  structure  of  FCSFIL.  Modifica¬ 
tions  were  also  made  to  the  output  formats  of  past  forecast  and  skill 
score  reviews. 

1.  Gerlach,  A.  M. ,  (ed) ,  1980a:  Computer-based  weather  research,  AFGL-TR- 

80-0069,  pp.  49-73. 

2.  Gerlach,  A.  M. ,  (ed) ,  1980b:  Development  or  automated  objective  meteor¬ 

ological  techniques.  AFGL-TR-81 -0017 ,  pp.  48-51. 


Fig.  1  shows  examples  of  the  format  used  in  listing  forecasts.  The 
primary  changes  involved  computation  of  vector  errors  for  wind  forecasts 
and  the  inclusion  of  the  variable  forecast  length  N.  Examples  of  skill 
score  reviews  for  the  first  two  forecast  lengths  are  presented  in  Fig.  2. 
Changes  for  Phase  II  included  elimination  of  three  types  of  skill  scores 
which  had  proved  uninformative  durinq  Phase  I  and  addition  of  two  new 
scores.  Percent  correct,  AWS  (Air  Weather  Service)  Srcill  score,  and  biar 
were  replaced  by  RMS  (root  mean  square)  error  and  the  Heidke  score.  P- 
score ,  cumulative  P-score,  and  mean  absolute  error  were  retained  from 
Phase  I . 

III.  Forecast  Procedure 

Durinq  April  -  May,  1981,  six  significant  weather  situations,  rang¬ 
ing  from  an  advancing  cloud  deck  to  developing  sea  breezes,  were  desig¬ 
nated  forecast  days.  Two  shifts  (morning,  afternoon),  of  three  fore¬ 
casters  each,  were  responsible  for  monitoring  and  preparing  hourly  fore¬ 
casts  for  each  situation.  At  their  disposal  were  the  entire  Me I DAS  now¬ 
casting  capability,  ha  1 f-hourly  one-mile  resolution  GOES  imagery,  the 
full  McIDAS  display  and  enhancement  capability,  OZ  and  ] 2Z  LFM  Model  Out¬ 
put  Statistics  (MOS)  forecasts,  and  OZ  and  12Z  LFM  products.  One-,  two-, 
and  four-hour  forecasts  of  total  cloud  amount  and  wind  speed  and  direction 
were  formulated  and  entered  for  Boston  and  a  second  station  in  the  north¬ 
east  U.S.  chosen  according  to  prevailing  conditions.  A  forecast  was  also 
prepared  and  entered  for  a  period  chosen  so  as  to  verify  on  the  nearest 
MOS  forecast  valid  time.  Forecasters  were  responsible  for  their  individ¬ 
ual  forecasts  but  group  discussions  were  encouraged  during  the  analysis 
period.  The  OZ  and  12Z  MOS  forecasts  and  forecasts  based  on  morning  and 
afternoon  climatology  for  Boston  were  entered  as  controls.  An  automated 
procedure  was  used  the  next  day  to  verify  the  forecasts  and  calculate 
skill  scores. 

Each  forecaster  selected  a  specific  facet  of  the  interactive  fore¬ 
cast  facility  for  concentrated  study  throughout  the  experiment.  Satel¬ 
lite  imagery  utilization,  new  data  sources,  improved  analyses  and  dis¬ 
plays,  and  forecast  aids  were  the  technical  areas  chosen,  based  on  indi¬ 
vidual  expertise  and  system  capabilities.  Advantages  and  disadvantages 
of  the  current  faciLity  were  examined,  and  potential  improvements  were 
considered. 
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IV.  Resulting  Improvements 


a.  Satellite  Imagery  Utilization 

One  area  of  interest  was  the  application  ui  real-time  geosta¬ 
tionary  satellite  imagery  to  mesoscalo  forecasting  ot  cloud  cover  and 
winds.  Loops  of  high  resolution  imaees  were  found  to  Le  invaluable  in 
predicting  the  motion  of  cloud  systems  and  in  timing  the  arrival  of  breaks 
and  cloud  bands  at  the  forecast  site.  Cumulus  growth  on  several  occasions 
signaled  the  development  of  a  sea  breeze  front  in  coastal  regions  where 
surface  reports  were  not  available.  However,  two  problems  with  the  satel¬ 
lite  display  capability  became  immediately  evident. 

The  more  serious  difficulty  was  inadequate  navigation,  or  earth 
location,  of  images.  Without  hiqhly  accurate  navigation  it  was  impossible, 
for  example,  to  judge  precisely  the  distance  between  the  forecast  scation 
and  an  advancing  cloud  edge.  Inadequate  navigation  also  resulted  in  appar¬ 
ent  motion  of  land  beneath  the  clouds.  Problems  were  encountered  in  at¬ 
tempting  to  separate  actual  cloud  motion  from  false  motions  induced  by 
navigation  error,  particularly  in  cases  where  no  land  features  were  visible 
in  the  images. 

To  improve  navigat. ional  accuracy,  copies  of  navigation  software 
(both  library  and  load  modules)  were  obtained  “'rom  the  apace  .Science  and 
Engineering  Center,  University  of  Wisconsin,  developers  of  the  McIDAS  sys¬ 
tem.  These  programs  were  modified  to  ntuko  them  compatible  with  the  AFGL 
McICAS  and  entered  into  tiie  system.  The  new  nav’ cut  ior.  technique  is  fully 
automated,  requiring  a  single  key-in  and  a  maximum  of  seconds  of  execu¬ 
tion  time.  This  replaces  a  process  which  required  approximately  four  hours. 
Imagery  loops  of  clearly  identifiable  land  features  remain  motionless  ever, 
at  full  resolution,  indicating  the  marked  improvement  in  navigational  ac¬ 
curacy. 

The  second  problem  noted  in  using  satellite  imagery  was  the  oc¬ 
casional  degradation  in  picture  contrast.  Programs  for  enhancing  the 
imagery  to  produce  greater  black  and  white  contrast  were  available,  but 
tended  to  be  time-consuming  and  awkward  to  use.  A  new  enhancement  pro¬ 
gram  was  written  which  scans  the  image,  determines  the  maximum  and  min¬ 
imum  brightness  values  present,  and  applies  a  contrast  stretch  that  opti¬ 
mizes  use  of  the  full  gray  scale  of  the  video  system.  The  dramatic  dif¬ 
ference  made  by  the  optimum  enhancement  program  is  illustrated  in  Figs.  3 


and  4.  Fig.  3  shows  an  unenhanced  infrared  image.  Fig.  4  shows  the  same 
image  after  enhancement.  Contrast  between  different  cloud  types  is  sharp¬ 
ened  as  well  as  the  cloud/land  and  land/ocean  interfaces.  The  enhanced 
image  even  reveals  variation  in  the  sea  surface  temperature  in  the  vicin¬ 
ity  of  the  Gulf  Stream.  The  speed  ai  1  simplicity  of  the  program  make  it 
a  very  useful  aid  to  the  forecaster  in  sharpening  details  of  cloud  fea¬ 
tures. 

b.  Increased  McIDAS  Data  Base 
(1)  FOUS  Bulletins 

Another  area  of  interest  during  MFE  Phase  II  was  evaluation  of 
the  current  McIDAS  non-satellite  meteorological  data  display.  One  goal 
of  this  evaluation  was  to  develop  a  series  of  routines  which  archive  and 
display  new  types  of  meteorological  data.  Initially,  only  surface  and 
upper  air  (RAOB)  reports  were  available  to  forecasters  using  McIDAS  for 
plotting  and  analysis  purposes.  Software  was  therefore  developed  to  ac¬ 
cess  and  display  new  non-satellite  data  types  available  through  land  base 
telephone  lines. 

To  provide  the  mesoscale  forecast  group  with  reliable  station- 
specific  forecast  guidance,  software  routines  were  written  which  decode, 
store  and  display  the  U’M-II  derived  United  States  Forecast  (FOUS)  bul¬ 
letins.  These  bulletins  are  available  via  the  McIDAS  conventional  data 
port  (weather  band  line  604)  . 

Three  series  of  FOUS  bulletins  are  available  for  display: 

F0US12  Model  Output  Statistics  (MOS) 

FOUS50-57  Trajectory  Model  Output 

FOUS60-78  Initial  LFM  Guidance  Output. 

Forecasts  contained  in  the  FOUS12  bulletin  are  highly  dependent 
on  output  from  the  LFM-II  model.  All  forecasts  are  produced  by  using  the 
MOS  technique  which  determines  a  statistical  relationship  between  the  pre- 
dictand  and  variables  from  the  numerical  model  at  some  prediction  time. 
Forecasts  being  saved  include  probability  of  precipitation,  categorical 
probability  of  quantitative  precipitation,  thunderstorm  probability,  con¬ 
ditional  probability  of  precipitation  type,  probability  of  heavy  snow, 
max\mum/mi nimum  temperature,  3-hour ly  surface  temperature  and  dew  point, 
surface  wind  speed  and  direction,  categorical  probabi 1 i ty  of  cloud  amount, 
categorical  probability  of  ceiling  height  and  visibility,  and  categorical 


probability  of  non-precipitating  obstructions  to  visibility. 

The  bulletins  designated  as  FOUS50-57  contain  24-hour  forecasts 
of  temperature,  dewpoint,  and  George's  K-index  for  each  station.  They 
also  contain  the  positions  (latitude,  longitude,  pressure  level),  at  6 
hour  intervals,  for  air  parcels  that  will  be  located  at  a  station  (sur¬ 
face,  350  mb,  and  700  mb)  at  the  end  of  the  24-hour  forecast  period. 

These  positions,  which  define  three  air  parcel  trajectories,  are  computed 
from  smoothed  wind  forecasts  generated  by  the  U-'M  model.  Knowledge  of 
air  parcel  movement  along  with  the  subsequent  evolution  of  the  temperature 
and  moisture  fields  provides  necessary  clues  to  the  development  of  con¬ 
vective  activity  and  precipitation. 

The  FOUS60-78  bulletins  contain  the  numerical  model  analysis 
and  fore casts  for  stations  in  the  U.  S.  and  adjacent  waters.  Forecasts 
are  provided  in  6-hourly  intervals  from  6-48  hours  after  the  model  run 
time  (either  0  or  12  GMT).  The  forecasts  are  relative  humidity,  vertical 
velocity  at  700  mb  Level,  lifted  index,  1000-500  mb  thickness,  direction 
and  speed  of  the  boundary  layer  wind,  mean  potential  boundary  layer  temper 
ature,  sea  level  pressure,  and  the  6-hourly  accumulated  precipitation. 
Figs.  5  through  7  illustrate  station  distribution  for  each  bulletin. 

The  FOUS  decoder  package  was  designed  to  be  as  transparent  to 
the  McIDAS  user  as  possible.  Much  of  the  software  is  activated  through 
internal  branching  from  the  McIDAS  data  line  handler  (SVCARX)  to  the  main 
FOUS  decoder  program  SVCFOS.  A  detailed  flow  chart  of  the  decoder  package 
is  given  in  Fig.  8. 

The  primary  function  of  SVCFOS  is  to  identify  the  test  received 
from  SVCARX  by  type  (either  F0US12,  50-57,  or  60-78)  and  transfer  process¬ 
ing  to  the  appropriate  subroutine.  Each  of  these  subroutines  contains 
code  specific  to  the  format  of  its  intended  bulletin.  It  is  designed  to 
sense  text  length  as  well  as  forecast  time(s)  and  file  storage  locations. 

Because  of  the  variety  in  formats  among  these  bulletins,  informa 
tion  decoded  from  the  test  is  stored  in  one  of  three  disk  data  files. 

Data  files  contain  station  3-letter  identifiers,  latitudes,  and  longitudes 
plus  a  tile  directory  and  blocked  storage  for  up  to  four  days  of  forecasts 
The  files  continually  update  themselves,  discarding  the  oldest  information 
to  provide  space  for  new  forecasts  as  they  are  decoded. 

The  display  program  FOSLIS  allows  users  access  to  information 
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stored  in  any  of  the  FOUS  data  files.  Listinqs  of  forecasts  for  any  data 
period  in  the  files  are  presented  in  the  NWS  format  on  the  CRT.  FC3LIS 
also  displays  file  directory  information  such  as  forecast  times  in  the 
file  and  the  number  of  stations  available  for  each  forecast  time.  Exam¬ 
ples  of  bulletin  display  formats  are  qiven  in  Fiqs.  9  and  10. 

(2)  Manually  Diqitized  Radar  (MDR)  Data 

SASC  meteoroloqists  modified  and  implemented  a  McIDAS  proqram 
set  partially  developed  by  the  University  of  Wisconsin,  which  decodes, 
stores,  and  displays  Manually  Digitized  Radar  tMDR)  code.  The  code,  which 
is  issued  hourly  by  NWS  radar  sites,  is  available  to  McIDAS  via  the  FAA 
hiyh  speed  data  link  (WB-604) . 

The  radar  code  is  made  up  of  two  major  sections:  an  SD  section 
which  gives  radar  data  in  an  azimuth-range  format  and  the  digital  section 
which  contains  echo  locations  and  intensities  in  digitized  form.  It  is 
the  digital  section  of  the  radar  code  which  is  processed  by  McIDAS. 

The  digital  section  is  based  upon  a  qrid  system.  Maximum  ob¬ 
served  Digital  Video  Integrator  and  Processor  (DVIP)  levels  are  reported 
for  each  grid  box  containing  echoes  of  moderate  or  greater  intensity. 

DVIP  levels  and  their  corresponding  intensities  are  listed  in  Table  1. 

Data  entered  in  the  digital  section  are  encoded  in  one  or  more 
qroups.  Each  group  starts  with  the  two-letter  identifier  of  the  grid 
box  in  the  row  containing  the  most  westerly  echo.  The  first  number  fol¬ 
lowing  the  box  identifier  is  the  DVIP  level  in  that  box.  Subsequent  num¬ 
bers  are  intensities  in  the  sequentially  adjacent  boxes  in  that  row. 

This  process  is  repeated  until  there  are  no  more  rows  containing  echoes 
to  be  encoded.  The  digital  section  is  limited  to  a  125  nm  ranqe.  Any 
echo  outside  this  limit  will  not  be  encoded. 

Fig.  11  illustrates  echoes  on  a  radar  scope  over  which  the  125 
nm  sweep  and  the  station  grid  system  have  been  laid.  The  complete  radar 
code  report  for  this  situation  is  listed  below  it.  The  start  of  the  dig¬ 
ital  section  is  marked  by  the  upward  turned  arrow.  Fiq.  12  illustrates 
the  encoded  DVIP  levels  as  they  would  appear  in  their  grid  box  positions. 

The  qrid  used  with  the  digital  code  is  a  subgrid  of  the  LFM-1I 
model  grid.  It  has  1/4  the  mesh  length  of  the  LFM  model  (boxes  are  ap¬ 
proximately  22  nm  on  a  side  and  have  an  area  of  484  sq  nm) .  Grid  boxes 
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TABLE  1.  RAINFALL  RATE  ECHO  INTENSITIES 
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Fiq.  11.  Example  of  Radar  Display  and  Radar  Code. 


are  located  with  respect  to  the  radar  site  such  that  box  MM  will  always 
contain  the  radar  antenna. 

Fig.  13  shows  the  digital  radar  code  grid  for  Worcester,  MA  (the 
skew  of  the  grid  is  due  to  the  curvature  of  the  earth) .  Table  2  lists 
the  National  Weather  Service  network  radars,  FAA/ARTC  radars,  U.S.  Air 
Force  radars,  and  NWS  local  warning  radars  from  which  McIDAS  receives 
digital  data. 

'r'R  software  is  desiqned  to  run  almost  independently  of  McIDAS 
users.  A  :  the  processing  routines  (decoding  and  store  programs)  are 
initiated  by  internal  cues  from  one  program  to  the  next.  Fig.  14  gives 
.i  detailed  flow  chart  of  the  MDR  decoder. 

Initialization  occurs  when  SVCARX  detects  an  ASCII  coded  string 
containing  MDR  information.  'Die  data  handler  sends  all  necessary  informa¬ 
tion  to  the  main  MDR  routine  SVCMDR,  which  searches  each  radar  report  for 
digitized  code.  If  found,  the  digitized  report  is  stored  within  the  MDR 
recovery  file  (MDRFI1.)  .  When  the  end  of  a  MDR  text  string  is  sensed,  the 
final  processing  program  MDRMAP  retrieves  the  code  stored  in  MDRFIL,  trans 
forms  the  digitized  field  to  earth  latitude  and  longitude  coordinates, 
and  calculates  the  McIDAS  MDR  display  qrid  point  field,  which  is  then 
stored  in  disk  qrid  file  GRIDOO. 

The  display  program  MOVMDR  allows  McIDAS  users  access  to  the  MDR 
grid  point  fields.  The  operator  simply  specifies  the  time  and  location  of 
the  data  to  be  displayed.  McIDAS  then  determines  map  specifications  and 
projections  for  the  specified  location  (map  specifications  include  region¬ 
al,  state,  city,  center  point,  and  satellite  overlay  options)  and  gener¬ 
ates  a  contoured  MDR  field  usinq  the  MDR  grid  point  data.  Figs.  15-17 
illustrate  the  ways  data  may  be  displayed  using  various  map  specifications 

c.  Improvements  In  Conventional  Data  Displays 

In  preparing  short-ranqe  terminal  forecasts,  regional  plots  of 
weather  parameters  are  valuable  aids.  Precise  tracking  of  dry  lines, 
wind  shifts,  and  precipitation  areas,  for  example,  can  keep  the  forecaster 
abreast  of  rapidly  changing  situations.  However,  the  density  of  data  plot 
ted  by  McIDAS  routines  IASVCA  and  WISPLT  was  significantly  less  than  the 
density  of  stations  reporting  data.  Figs.  18a  and  19a  show  that  this  re¬ 
duction  was  done  to  prevent  overlapping  of  reports  on  the  mat'  projections 


TABLE  2.  SOURCES  OF  DIGITAL  RADAR  DATA  FOR  MCIDAS 
NATIONAL  WEATHER  SERVICE  NETWORK  RADARS 


Alliance,  Nebraska 
Amarillo,  Texas 
Apalachicola,  Florida 
Athens,  Georgia 
Atlantic  City,  New  Jersey 
Binghamton,  New  York 
Bristol ,  Tennessee 
Brownsville,  Texas 
Brunswick,  Maine 
Buffalo,  New  York 
Cape  Hatteras,  North  Carolina 
Centreville,  Alabama 
Charleston,  South  Carolina 
Chatham,  Massachusetts 
Cincinnati,  Ohio 
Daytona  Beach,  Florida 
Des  Moines,  Iowa 
Detroit,  Michigan 
Evansville,  Indiana 
Fargo,  North  Dakota 
Galveston,  Texas 
Garden  City,  Kansas 
Grand  Island,  Nebraska 
Hondo,  Texas 
Huron,  South  Dakota 
Jackson,  Mississippi 
Kansas  City,  Missouri 
Key  West,  Florida 


FAA/ARTC 

Albuquerque,  New  Mexico 
Amarillo,  Texas 
Ashton,  Idaho 
Battle  Mountain,  Nevada 
Boise,  Idaho 
Boron,  California 
Cedar  City,  Utah 
El  Paso,  Texas 
Klamath  Falls,  Oregon 
Las  Vegas,  Nevada 
Lovell,  Wyoming 


Lake  Charles,  Louisiana 
Limon,  Colorado 
Little  Rock,  Arkansas 
Longview,  Texas 
Marseilles,  Illinois 
Medford,  Oregon 
Memphis,  Tennessee 
Miami,  Florida 
Midland,  Texas 
Minneapolis,  Minnesota 
Missoula,  Montana 
Monett,  Missouri 
Nashville,  Tennessee 
Neenah,  Wisconsin 
New  York  City,  New  York 
Oklahoma  City,  Oklahoma 
Patuxent  River,  Maryland 
Pensacola,  Florida 
Pittsburgh,  Pennsylvania 
Sacramento,  California 
Slidell,  Louisiana 
Stephenville ,  Texas 
St.  Louis,  Missouri 
Tampa,  Florida 
Volens,  Virginia 
Waycross,  Georgia 
Wichita,  Kansas 
Wilmington,  North  Carolina 


RADARS 

Mesa  Rica,  New  Mexico 
Mt .  Laguna,  California 
Paso  Robles,  California 
Phoenix,  Arizona 
Rock  Springs,  Wyoming 
Salem,  Oregon 
Salt  Lake  City,  Utah 
San  Pedro,  California 
Seattle,  Washington 
Silver  City,  New  Mexico 
Spokane,  Washington 


U.S.  AIR  FORCE  RADARS 


Laughlin  AFB,  Texas 


Loring  AFB,  Maine 
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TABLE  2  (cont.) 


NATIONAL  WEATHER  SERVICE  LOCAL  WARNING  RADARS 


Abilene,  Texas 
Akron,  Ohio 
Albany,  New  York 
Alpena,  Michigan 
Atlanta,  Georgia 
Augusta,  Georgia 
Austin,  Texas 
Baton  Rouge,  Louisiana 
Beckley,  West  Virginia 
Billings,  Montana 
Bismarck,  North  Dakota 
Burlington,  Vermont 
Charleston,  West  Virginia 
Charlotte,  North  Carolina 
Cheyenne,  Wyoming 
Cleveland,  Ohio 
Columbia,  Missouri 
Columbia,  South  Carolina 
Columbus,  Georgia 
Columbus,  Ohio 
Concordia,  Kansas 
Corpus  Christ i,  Texas 
Duluth,  Minnesota 
Erie,  Pennsylvania 
Fort  Smith,  Arkansas 
Fort  Wayne,  Indiana 
Goodland,  Kansas 
Harrisburg,  Pennsylvania 
Hartford,  Connecticut 
Houghton  Lake,  Michigan 
Huntsville,  Alabama 
Indianapolis,  Indiana 
Jackson,  Kentucky 
Las  Vegas,  Nevada 


Los  Angeles,  California 
Louisville,  Kentucky 
Lubbock,  Texas 
Macon,  Georgia 
Madison,  Wisconsin 
Marquette,  Michigan 
Meridian,  Mississippi 
Mobile,  Alabama 
Moline,  Illinois 
Montgomery.  Alabama 
Muskegon,  Michigan 
Norfolk,  Nebraska 
North  Platte,  Nebraska 
Omaha,  Nebraska 
Phoenix,  Arizona 
Portland,  Oregon 
Raleigh,  North  Carolina 
Rapid  City,  South  Dakota 
Rochester,  Minnesota 
San  Angelo,  Texas 
San  Juan,  Puerto  Rico 
Shreveport,  Louisiana 
Sioux  Falls,  South  Dakota 
Springfield,  Illinois 
Topeka,  Kansas 
Tulsa,  Oklahoma 
Victoria,  Texas 
Waco,  Texas 
Waterloo,  Iowa 
West  Palm  Beach,  Florida 
Wichita  Falls,  Texas 
Williston,  North  Dakota 
Worcester,  Massachusetts 


S  V  C  A  R  X 

MclDAS  conventional  data  handier 

SVCMDR 

MDR  pre-processing  program 

^7 

MDRFIL 

storage  of  digitized  radar  code 

MDRMAP 

conversion  from  MDR  code 
to  MclDAS  grid  point  field 

^7 

G  R  I  D  0  0 

MclDAS  grid  point  file 

^7 

M  0  V  M  D  R 

MDR  display  program 


Fig.  14.  MDR  Decoder  Flow  Chart. 


HDR  SFC  131700 


Fig.  15.  MDR  Regional  Display  for  Eastern 
United  States. 


provided.  In  order  to  plot  all  available  data,  a  smaller  map  projection 
was  devised,  to  be  centered  on  any  of  the  1103  stations  in  the  McIDAS 
surface  data  file.  Figs.  18b  and  19b  illustrate  the  resultant  increase 
in  plotted  data,  approximately  threefold,  which  allows  more  precise  lo¬ 
cation  of  discontinuities. 

The  capability  to  display  four  analyses,  one  in  each  quadrant 
of  the  CRT  monitor,  was  recognized  as  an  important  feature  of  interactive 
forecast  sy -terns.  Four-panel  analysis  of  various  parameters,  such  as 
dat"  from  a'ious  levels  in  the  atmosphere  or  data  from  various  time  per¬ 
iods,  facilitates  rapid  visual  integration  of  pertinent  data.  To  provide 
data  pi  'tting  overlay  compatibility  as  well  as  four-panel  plotting  capa¬ 
bility,  the  graphics  plotting  routines  IASVCA  and  IARAOB  were  rewritten 
to  include  that  option  (see  Fig.  20). 

When  these  plotting  capabilities  were  used  operationally  it  was 
noted  that  stations  with  incomplete  reports,  i.e.,  temperature,  dew  point 
temperature,  or  pressure  reported  as  missing,  would  not  be  plotted  at  all. 
The  missing  data  checks  in  the  plotting  routines  were  revised  to  permit 
plotting  of  reported  data  in  such  incomplete  reports.  The  data  listing 
routine  SVCADI  was  also  revised  to  list  incomplete  station  reports  which 
were  previously  left  blank  (see  Fig.  21)  . 

d.  Interactive  Forecast  Aids 

In  response  to  the  need  for  forecast  aids  noted  during  the  ex¬ 
periment,  two  menu-driven  quidance  routines  were  added  to  the  MFF :  a 
cloud  band  movement  decision  tree  algorithm  developed  at  Pennsylvania 
State  University^  and  an  automated  version  of  the  severe  weather  check¬ 
list  developed  by  the  National  Weather  Service  Forecast  Office  at  Portland, 
4 

Maine.  Both  routines  ask  a  series  of  yes/no  questions  designed  to  guide 
the  forecaster  through  an  analysis  of  the  current  weather  situation.  The 
cloud  band  decision  tree  helps  isolate  the  key  factor  governing  the  mo¬ 
tion  of  the  cloud  band,  i.e.,  confluence  zones  or  vorticity  patterns. 

3.  Cahir,  J.,  J.  Norman  and  D.  Lowry,  1981:  Use  of  a  real-time  computer 

graphics  system  in  analysis  and  forecasting.  Monthly  Weather  Review, 

109,  486-488.  "  . 

4.  Gulezian,  D. ,  1980:  Severe  weather  checklist  used  at  National  Weather 

Service  Forecast  Office,  Portland,  Maine.  Bulletin  of  the  American 

Meteorological  Society,  61 ,  1592-1599. 


4 

1 


IASVCA  Regional  Plot  of  Temperature. 


Ml 


Fig.  19a.  WISPLT  Regional  Plot  of  Temperature. 


Fig.  19b.  WISPLT  Detailed  Regional  Plot  of  Temperature. 


Four-Panel  Plot  of  Current  Weather 


When  this  key  factor  is  determined,  a  brief  message  advising  forecast  pro¬ 
cedure  is  displayed.  The  checklist  tabulates  yes  answers  provided  by  the 
forecaster  and  predicts  whether  convective  weather  will  occur  and  indi¬ 
cates  whether  or  not  it  will  reach  severe  levels.  both  routines  are  areu 
forecasting  aids  that  can  be  applied  to  a  particular  station, 
e.  Automated  Analysis  and  Forecast  Routines 

When  a  previously  calculated  analysis  was  redisplayed  on  the 
WRRRM  graphics,  AFGL  surface  and  upper  air  analyzing  routines  SACf'KT  and 
RPCONX  had  to  collect  the  raw  data,  re-int:erpolate  them  to  grid  ■ 
and  then  contour  those  re-gridded  data.  To  speed  up  this  redisplay  \  re¬ 
cess,  grid  files  were  added  to  the  McIDAS  data  base  to  store  interpolated 
data.  Software  developed  at  the  University  of  Wisconsin  to  write  gridded 
data  into  these  files  and  to  display  saved  analyses  was  modified  and  added 
to  the  nowcasting  capabil ity . 

With  gridded  data  saved  in  this  manner,  processing  and  use  of 
the  grids  not  previously  possible  could  now  be  accomplished.  A  bilinear 
interpolation  technique  to  return  the  field  value  at  any  point  in  the  an¬ 
alysis  was  developed  and  added  to  PLTMAP.  This  option,  when  used  with  a 
two-dimensional  approximation  of  the  LFM-derived  Trajectory  Model  Fore¬ 
cast,  provided  a  means  to  forecast  any  conservative  property  of  the  at¬ 
mosphere  at  the  surface,  850,  and  700  millibar  levels  for  12-  and  24-hour 
time  periods. 

Program  GRIDWK,  developed  at  the  University  of  Wisconsin,  was 
modified  and  added  to  McIDAS  to  manage  and  manipulate  the  grids  within 
each  file.  Options  to  change  extreme  values  within  the  grid  permit  crude 
data  editing  and  changes  to  the  analysis  without  affecting  the  original 
raw  data.  The  ability  to  transfer  grids  from  the  file  of  one  terminal  to 
another  provides  common  access  to  any  analysis.  In  response  to  a  need 
to  archive  MDR  data,  which  must  be  transferred  to  a  terminal’s  grid  file 
to  be  displayed,  save  and  restore  options  were  developed  for  GRIDWK.  Up 
to  28  eastern  U.S.  0.5°  resolution  grids  of  MDR  data  can  be  saved  or  re¬ 
stored  by  these  routines. 

Access  to  gridded  data  also  allowed  the  combination  of  surface 
and  upper  air  data  in  a  more  accurate  1000-500  millibar  thickness  analysis. 
Since  a  1000  mb  level  report  is  not  received  from  a  RAOB  station  if  its 
uncorrected  surface  pressure  is  less  than  1000  mb,  previous  analysis 
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routines  which  were  limited  to  upper  air  data  only  contained  large  re¬ 
gions  with  missing  data.  By  using  a  OZ  or  12Z  surface  pressure  analysis 
and  converting  it  to  a  1000  mb  height  grid,  a  quick  and  reasonable  ap¬ 
proximation  to  the  thickness  field  can  be  calculated. 

Writing  gridded  data  to  a  file  permitted  the  addition  to  McIDAS  of 
two  routines  developed  at  the  University  of  Wisconsin  to  aid  the  short- 
range  forecaster.  PRERAB,  a  trajectory  model  forecast  routine,  and 
SETRAB,  a  vertically  averaging  wind  field  analysis  routine,  were  modi¬ 
fied  for  compatibility  with  the  AFGL  McIDAS.  A  vertically-averaged, 
density -weighted  wind  field  using  any  or  all  mandatory  level  data  is  cal¬ 
culated  to  produce  u  and  v  component  wind  field  grids.  Using  these  grids 
and  an  initial  field  grid,  point  or  station  forecasts  and  area  forecasts 
based  on  either  advection  or  trajectory  methods  can  be  computed  for  any 
time  period. 

V.  Conclusion 

The  speed  and  flexibility  of  McIDAS  in  analyzing  and  displaying  the 
data  base  available  have  only  begun  to  be  exploited.  The  addition  of  LFM 
forecast  products,  simple  forecast  techniques,  and  menu-driven  guidance 
programs  allows  the  terminal  forecaster  to  assess  the  regional  impact  of 
weather  systems.  Expansion  of  the  McIDAS  data  base  to  include  MDR  data 
and  improvements  in  GOES  imagery  navigation  provide  important  data  to 
the  mesoscale  forecaster  for  modifying  this  regional  picture  and  tailor¬ 
ing  it  to  his  local  area.  While  modification  of  McIDAS  may  not  yield 
the  ideal  terminal  forecast  system,  the  changes  made  do  point  out  the 
direction  in  which  the  development  of  such  a  system  should  proceed. 

B.  AFGL  WEATHER  TEST  FACILITY  SUPPORT 

I.  Introduction 

The  Weather  Test  Facility  (WTF)  at  Otis  AFB  of  the  Air  Force  Geo¬ 
physics  Laboratory  includes  towers  equipped  with  weather  instruments  at 
selected  heights  which,  along  with  ground  based  instruments,  sample  wind, 
temperature ,  visibility,  precipitation,  cloud  base  height  and  other  weath¬ 
er  parameters. 

To  facilitate  data  collection,  the  Modular  Automated  Weather  System 
(MAWS)  was  installed  at  Otis  to  accept  signals  from  the  various  sensors, 
create  one-minute  mean  digital  values  and  send  these  values  over  tele- 
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phone  long  lines  to  AFGL  where  they  are  recorded  onto  magnetic  tape. 


II.  Data  for  Sensor  Evaluation 

Under  the  contract  SASC  was  involved  in  computerized  data  process¬ 
ing  and  evaluation  of  new  or  improved  sensors  from  MAWS  data  tapes.  The 
sensors  included  Rosemount  wind  sets  and  such  visibility  meters  as  a  mod¬ 
ified  Fumosens,  Wright  forward  scatter  visibility  meter,  an  AFGL-devel- 
oped  side- looking  forward  scatter  meter  referred  to  as  a  "lookaway , "  and 
a  Wriqht  Fog  15  visibility  sensor. 

The  processing  was  accomplished  by  various  versions  of  a  CDC 
ter  program,  MAWSE,  which  displays  outputs  of  the  sensors  in  convenient 
form  for  analysis.  Statistical  correlation  between  the  data  from  the 
instrument  being  tested  and  those  from  one  of  known  reliability  was  com¬ 
puted  and  printed  along  with  simple  statistical  measures  of  the  two. 
Additional  statistical  information  was  displayed  on  graphical  represen¬ 
tations  used  to  high, light  similarities  and  differences  between  instru¬ 
ments  . 

Another  effort  in  sensor  evaluation  involved  examining  periods  of 
significant  wind  shear.  A  wind  shear  event  was  defined  as  a  period  of 
time  when  tire  wind  shear  between  the  200  ft  level  and  the  ground  was 
greater  than  four  meters  per  second  for  a  period  of  more  than  four  min¬ 
utes'  duration.  Episodes  of  wind  shear  were  isolated  from  the  minute  to 
minute  wind  data.  Tire  shear  magnitude  of  a  given  event  was  examined  in 
relation  to  the  surface  wind  during  that  event  and  in  relation  to  the 
temperature  differences  occurring  between  the  two  heights  during  the  event. 

III.  Verification  of  Prediction  Techniques 

A  further  Air  Force  concern  was  short-term  prediction  techniques  for 
weather  parameters  critical  to  aviation.  Computer  program  VERIFY,  designed 
to  be  the  final  software  module  in  AFGL's  low  cloud  observation  and  pre¬ 
diction  system,  was  completed.  The  program  performed  a  thorough  analysis 
of  all  probability  and  cateqorical  forecasts  based  on  cloud  ceiling  data 
retrieved  from  three  rotating  beam  ceilometers  (RBC)  and  wind  and  vis¬ 
ibility  data  as  processed  by  the  MAWSE  program.  Forty  episodes  durinq 
which  clouds  were  reported  below  6000  ft  for  at  least  three  hours  were 
selected.  Approximately  20,000  one-minute  observations  were  processed. 


Verification  statistics  were  calculated  for  15  prediction  techniques 
which  fell  into  two  groups.  The  first  group  consisted  of  probability 
forecasts ,  ba-e.i  on  KEEP  equations.  REK!  (Regression  Estimation  of  Event 
;  rob.ibi  1  it  ;es'  .  ab  ulates  probabilities  of  being  within  certain  categor¬ 
ies.  The  model  employed  placed  cloud  data  into  five  categories:  <200  ft, 
-jo-500  ft,  "0v'-lv00  ft,  1000-3000  ft,  and  >3000  ft  or  no  ceiling;  or 
into  two  categories,  <1000  ft  or  -  1000  ft.  Exceedance  probabilities 
were  prod  •  •:>.}  base-1  on  three  sets  of  KEEP  equations:  KEEP'S  using  only 
cloud  ob  vv.it ions,  REHP.Ol  using  the  predictor  set  in  which  the  last 
;  radio tor  selected  reduced  the  total  variance  by  at  least  one  percent, 
and  REEPll  where  wind  and  visibility  data  were  added  to  the  cloud  obser¬ 
vations.  Also,  an  equivalent  Markov  technique  and  persistence  were  used. 

The  second  group-  included  thresholding  techniques  such  as  maximum 
probability,  iterative,  using  REEP,  and  unit  bias,  climatology,  and  per¬ 
sistence  . 

All  the  exceedance  probabilities  and  categorical  forecasts  generated 
were  supplied  to  SAKC  to  be  used  in  program  VERITY  for  four  forecast 
lengths:  30,  GO,  120  and  180  minutes.  Statistical  tests  applied  to 

these  data  included  bias,  threat  score,  percent  correct,  Brier  score,  t- 
stutistic,  and  chi-square.  For  each  successive  ejisode  of  the  40  se¬ 
lected,  these  statistics  were  compiled  and  stored  in  a  random  access 
data  file  (RADF).  The  RAUF  retained  data  for  each  of  the  15  forecast 
techniques,  stratified  by  forecast  length,  episode  number,  and  eventually 
for  all  o(  isodes  combined.  The  KADI'  required  less  storage  than  a  sequen¬ 
tial  file,  and  provided  faster  and  more  efficient  retrieval  for  several 
persons  desir in  i  access  to  the  information,  whether  for  analysis,  or 
print,  or  plot  display. 

VERIFY  provided  a  summary  listing  of  all  statistics  calculated  and 
became  the  critical  module  whicn  provided  the  initial  forecast  process¬ 
ing  and  subsequent  display  of  verification  data. 

IV.  Fog  Model  Display 

A  program  was  written  to  produce  a  three-dimensional  display  of  the 
temporal  variability  of  the  drop  size  distribution  during  the  life  cycle 
<■>•*  i  tog  episode.  As  part  of  an  AFC.  I.  effort  to  dove  lot  improved  numer¬ 
ical  techniques  for  the  short-range  pt  edict  ion  of  tog,  the  drop-  size 
distribut  ion  data  were  obtained  at  the  WTF  with  a  r.MF  forward  scatter 


spectrometer  system.  a  graphics  display  package  developed  Ly  the  National 
Center  tor  Atmosi  he'  ic  Research  (NCAR)  was  extensively  revised  to  achieve 
the*  |  roi'Cr  angles  tor  displayin'.'  the  data  and  » 1  tl  <  - !  the  axes.  With 
add  it  ional  .outir.es,  displays  were  produce  i  wh  i  ch  i.  •  t.  •  r  t  •.«;  varia¬ 

tion  1  .  '.ho  droi  let  concent ;  at  ion  v  N )  with  dro:  ..i::o  md  time. 

: '  i  :.  !»:'  show.-  *  w’o  djs.pl ays  of  data,  .-ol  looted  it  n  and  a  x  above 
ground ,  during  i  f  ivc— hour  fog  episode  at  otis  Ai’i-. 
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Drop-Size  Distributions  From  Spectrometer  Probes 


SECTION  3. 

SATELLITE  METEOROLOGY 


A.  CLUSTERING  TECHNIQUE  FOR  SATELLITE  IMAGERY  ANALYSIS 


As  one  step  in  developing  an  automated  satellite  imagery  analysis 
scheme,  AFGL  developed  a  technique  using  infrared  radiative  count  values 
for  analyzing  data  samples  into  sub-regions  representing  various  cloud 
types. 1  Manual  classification  was  performed,  assigning  to  each  sub- 
region,  or  cluster,  a  cloud  type  based  on  the  scheme  in  Table  1. 

Texture  and  gradient  calculations  were  performed  on  the  data  in  each 
cluster.  It  was  expected  that,  with  a  sufficiently  large  number  of  sam¬ 
ples  from  each  of  the  eighteen  cloud  categories,  distinctive  attributes 
or  combinations  of  attributes  could  be  isolated  by  which  each  cloud  type 
could  be  identified  in  a  fully  automated  system.  All  data  selection  and 
analysis  were  performed  on  the  AFGL  C DC  6G00  computer  using  recorded 
satellite  imagery  from  the  ODIS  (Offline  Data  Ingest  System)  archive. 

TABLE  1.  CLOUD  CLASSIFICATION  SYSTEM. 


Cirrus 

Cumulus 

Othei 

1. 

Thin 

8. 

Sma 1 1  scattered 

12. 

Stratus 

2 . 

Moderate 

9. 

Sma  1 1 

13. 

Stratocumul us 

3. 

Thick 

10. 

Moderate 

14. 

Al tocumulus 

4. 

Thin  over  Cu 

11 . 

Large 

15. 

Al tostratus 

5. 

Moderate  over 

Cu 

It. 

Nimbostratus 

G. 

Thin  over  St 

1  7. 

Cumulonimbus 

7. 

Moderate  over 

St 

18. 

Clear 

SASC  developed  an  interactive  -adaptation  of  the  clustering  algorithm, 
called  CLUSTR,  for  use  on  McIDAS.  Documentation  of  the  command  options  is 
provided  in  Appendix  A  to  this  section.  The  interactive  approach  was  se¬ 
lected  to  facilitate  quick  and  easy  sample  selection.  Analysis  and  dis¬ 
play  options  were  later  incorporated,  followed  by  file  editing  and  manage¬ 
ment  capabilities. 

A  stepwise  procedure  was  established  for  data  collection.  Visible 
and  infrared  satellite  images  were  loaded  on  television  frames  and  exam¬ 
ined  for  suitable  cloud  formations.  After  a  sample  was  visually  selected, 
the  box  cursor  under  joystick  control  was  positioned  over  it,  and  the 
command  CL  was  entered.  This  command  initiated  CLUSTR,  which  extracted 

1.  Hawkins,  R.  S. ,  1980:  A  clustering  technique  for  Satellite  imagery 
analysis.  Preprints ,  8th  Conference  on  Weather  Forecasting  and 
Analysis,  AMS,  Boston;  pp.  135-118. 
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the  visible  and  infrared  data  arrays  corresponding  to  the  region  within 
the  cursor,  executed  the  clustering  algorithm,  computed  histograms  of 
texture  and  gradient  at  each  of  three  spatial  scales,  and  entered  the 
results  in  a  disk  file  called  CLSFIL.  Documentation  of  the  structure  of 
CLSFIL  is  provided  in  Appendix  B. 

A  message  was  sent  to  the  operator's  console  indicating  where  the 
sample  was  stored  in  the  file.  CLSFIL  has  a  maximum  capacity  of  ten 
sample  entries.  Samples  were  displayed  on  the  television  monitor  by  the 
LOAD  option  of  CLUSTR.  Fig.  1  is  an  example  of  the  type  of  display  re¬ 
sulting  from  this  command.  The  visible  imagery  sample,  blown  up  by  a 
factor  of  three,  appears  in  the  upper  left  quadrant  of  the  picture.  The 
corresponding  infrared  sample  is  shown  in  the  lower  left.  In  the  upper 
right  quadrant  is  a  representation  of  the  results  of  clustering.  Each 
of  the  four  clusters  derived  from  the  infrared  data  is  depicted  in  a 
different  gray  shade. 

Based  on  the  information  provided  in  displays  like  Fiq.  1,  manual 
cloud  classification  was  performed  and  entered  in  the  TYPE  option  of 
CLUSTR.  Texture  and  gradient  histograms  were  drawn  using  the  option  HIS. 
Each  curve  could  be  displayed  individually  or  in  combination  with  up  to 
five  others.  Fig.  2  shows  an  example  of  histograms  of  the  three  scales 
of  texture  for  the  third  cluster  of  an  imagery  sample.  On  the  McIDAS 
graphics  display  the  TEX  1  curve  and  annotation  appeared  in  red,  TEX  2  in 
green,  and  TEX  3  in  yellow.  Additional  curves  were  added  in  dashed 
lines  up  to  a  maximum  of  six  curves  per  display,  as  in  Fig.  3.  Labels, 
not  present  on  McIDAS,  have  been  added  in  Figs.  2  and  3  to  identify  each 
curve. 

Summaries  of  data  collected  for  individual  samples  or  for  all  entries 
in  the  file  were  available  through  the  SUM  option  of  CLUSTR.  Listings 
were  directed  either  to  the  operator's  console  or  to  the  line  printer  for 
hard  copies,  like  the  example  in  Fig.  4.  The  entry  number  in  these  list¬ 
ings  refers  to  position  within  CLSFIL,  and  the  sample  number  refers  to 
location  on  a  particular  visible/infrared  picture  pair  identified  by 
SSYYDDD  (satellite-year-date  code)  and  HHMMSS  (hour-minute--second  of  pic¬ 
ture  start  time) .  Cut  levels  and  Mean  IR  and  VIS  are  given  in  radiative 
count  values.  Cloud  type  assignments  correspond  to  the  numbers  given  in 
Table  1. 
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Sample  Display  Produced  by  CT.USTR. 


Fig.  2.  Histogram  Results  for  Three  Scales  of  Texture. 


ae  at  at 

UJ  HI  UJ 

aaa 

UULJ 


and  Gradient  Calculations. 


*»* 


entry  1 


ENTRY  2 


Entky  j 


EMI  WY  (4 


EMTHY  S 


AEGL  MCIOAS  ***  PROGRAM j  CLU3TR 


contents  op  cluster  file 

sample  2  SSYYD00s?881252  HHMMSS*! 30000 


CLUSTER 

I 

2 

3  4 

cut  Ltvi'L 

115  155 

205 

MO  OF  POINTS 

38 

42 

119  57 

MEAN  iw 

105.0 

131.6 

188.6  211.1 

mean  vis 

73.4 

100.9 

149.0  179.9 

CLOUD  TYPE 

18 

1 

4  17 

SAMPLE  1 

SSYY0[)0s288l254 

HHMMSS*130000 

CLUSTER 

1 

2 

3 

CUT  LEVEL 

105  115 

vO  of  points 

78 

51 

129 

MEAN  IR 

95.6 

109.2 

124.0 

MEAN  VIS 

82.3 

107.9 

152.3 

CLOUO  type 

18 

1  3 

14 

sample  2  SSYYD0DS2881254 

HHMM3S*1 30000 

CLUSTER 

1 

2 

CUT  LEVEL 

95 

VO  OF  POINTS 

2$  9 

17 

MEAN  IR 

80.9 

101.5 

MEAN  VIS 

66.5 

96.0 

CLOUD  TYPE 

16 

15 

SAMPLE  5  SSYYODDs 

2o81254 

MHMMSS=130000 

CLUSTER 

1 

2 

CUT  LEVEL 

95 

MO  OF  POINTS 

38 

216 

MEAN  IH 

92.8 

104.4 

MEAN  VIS 

69.3 

123.1 

CLOUO  TYPE 

18 

13 

Sample  4  SSYYQODb 

2881254 

MHMMSS*130000 

; 

CLUSTtR 

1 

vO  OF  POINTS 

256 

MEAN  IR 

91.4 

MEAN  VIS 

96.8 

CLOUD  TYPE 

12 

Fig.  4.  Summary  Listing  of  Contents  of  CLSFIL. 
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The  final  CLUSTR  display  option,  BOX,  was  used  to  identify  all  sam¬ 
ples  selected  from  an  image.  Figs.  5  and  6  show  boxes  over  visible  and 
infrared  satellite  images  respectively  with  labels  corresponding  to  sam¬ 
ple  numbers  in  entries  2-5  of  Fig.  4. 

When  CLSFIL  was  full  or  sample  selection  was  completed,  file  con¬ 
tents  were  saved  on  magnetic  tape  using  the  SAVE  option  of  CLUSTR.  The 
converse  function,  that  of  restoring  data  from  tape  to  disk  file,  was 
performed  by  the  RES  option.  File  editing  functions  consisted  of  INI, 
which  initialized  or  cleared  the  contents  of  CLSFIL,  and  DEL,  which  was 
used  to  delete  any  single  entry  from  the  file. 

B.  BRIGHTNESS  VARIATION  ANALYSIS 

In  an  effort  to  utilize  geostationary  satellite  data  for  determining 

cloud  cover  amounts,  SASC  scientists  developed  software  to  reduce  and  an- 

2 

alyze  archived  GOES  and  SMS  imagery  (Gerlach,  1980a)  .  The  conclusions 
reached  in  the  first  phase  of  the  effort  —  calculation  of  curves  relat¬ 
ing  cloud  cover  and  average  reflectivity  (Method  1)  —  were  reported  in 
Keegan  and  Niedzielski  (1981) 

This  follow-on  study  entailed  the  testing  of  these  relationships  on 
the  dependent  data  sample  (1977-1978  data)  to  assess  their  effectiveness. 

Results  showed  an  overall  67  percent  correct  classification  rate 
(see  Table  2) .  Weather  station  reports  of  cloud  amounts  were  used  as 
cloudtruth.  Detailed  study  showed  this  rate  to  be  heavily  weighted  by 
the  clear  and  overcast  categories,  which  account  for  31  percent  and  24 
percent  of  the  sample,  respectively.  Both  categories  were  classified 
correctly  nearly  80  percent  of  the  time.  Categories  of  scattered  and 
broken  cloud  cover  had  correct  classification  rates  ranging  from  20  per¬ 
cent  for  three-tenths  coverage  to  57  percent  for  seven-tenths  coverage. 
Testing  with  an  independent  data  sample  (1979  data)  yielded  essentially 
the  same  results. 

2.  Gerlach,  A.  M.  (ed) ,  1980a:  Computer-based  weather  research.  AFGL- 

TR-80-0069 ,  pp.  74-81. 

3.  Keegan,  T.  J. ,  and  M.  E.  Niedzielski,  1981:  Specification  of  cloud 

amount  over  local  areas  from  GOES  visual  imagery.  Environmental 
Research  Papers  No.  743,  AFGL-TR-81-0153 . 


TABLE  2.  COMPARISON  OF  CLOUD  AMOUNT  TECHNIQUES 


METHOD 

1 

2 

3 

4 

All  Cloud  Amounts 

67 

64 

75 

77 

Clear 

80 

79 

85 

92 

Overcast 

80 

78 

91 

91 

Correct  Classification  Rates  in  Percent 

Method  1  is  Average  Reflectivity,-  Method  2  is  Average 
Reflectivity,  Texture,  and  Maximum  Reflectivity;  Meth¬ 
od  3  is  Simple  Threshold  (fixed  threshold) ;  and  Meth¬ 
od  4  is  Variable  Threshold  (decision  tree  algorithm) . 

Two  important  conclusions  were  reached  as  a  result  of  these  tests. 
First,  average  reflectivity  could  generally  distinguish  cloud-free  and 
cloud-filled  fields  of  view  but  was  a  poor  parameter  for  determining 
scattered  and  broken  cloud  amounts.  Second,  the  data  remained  stable 
enough  from  several  satellites  (GOES  1,  GOES  2,  SMS  2)  to  yield  nearly 
constant  results. 

The  addition  of  other  reflectivity  characteristics  for  multivariate 
analysis  was  the  next  phase  of  the  study.  Relationships  between  cloud 
amount  and  average  reflectivity,  sample  texture  (standard  deviation),  and 
maximum  reflectivity  were  calculated  (Method  2)  and  the  testing  repeated 
on  both  dei'endent  and  independent  data.  Overall  results  of  classifica¬ 
tion  rates  decreased  slightly  to  64  percent  correct.  This  led  to  the 
subjective  conclusion  that  the  variability  of  reflectivity  characteris¬ 
tics  v.'ithin  each  cloud  cover  category  was  too  large  to  allow  regression 
techniques  to  define  the  relationships  clearly. 

As  a  result  of  the  indifferent  ability  of  regression  curves  to  spec¬ 
ify  cloud  cover,  a  different  method,  a  threshold  technique,  was  implement 
ed  (Method  3).  This  approach  is  similar  to  the  Air  Force  Global  Weather 
Central  cloud  analysis  program  30NEPH,  which  uses  location-dependent 

thresholds  for  cloud  detection  in  visible  data  on  polar-orbiting  satel- 

4 

lites. 

4.  Fye,  F.  K. ,  1978:  The  AFGWC  automated  cloud  analysis  model.  AFGWC 
Technical  Memorandum  78-002,  pp.  26-31. 


Using  information  deduced  from  analysis  of  all  clear  and  overcast 
cases,  a  cloud/no  cloud  threshold  reflectivity  of  18  percent  was  select¬ 
ed.  A  simple  addition  was  made  to  program  CENTER  to  count  the  number 
of  pixels  within  the  9x9  pixel  analysis  box  that  were  above  the  thresh¬ 
old  value  and  to  compute  a  percent  cloud  cover.  This  straightforward 
technique  resulted  in  a  75  percent  correct  classification  rate  of  cases 
with  no  tenuous  cirrus  clouds.  Clear  and  overcast  conditions  were  clas¬ 
sified  correctly  85  percent  and  91  percent  of  the  time,  respectively. 

Closer  analysis  of  misclassif ied  clear  cases  showed  that  sample  texture 
could  help  distinguish  cases  with  clouds  whose  areal  extent  was  somewhat 
below  image  resolution.  If  the  sample  texture  exceeded  2.5  percent,  a 
cloud/no  cloud  threshold  of  16.0  percent  was  indicated.  Incorporating 
this  information  into  a  simple  decision  tree  algorithm  in  CENTER  (Meth¬ 
od  4)  resulted  in  a  92  percent  correct  classification  rate  for  clear  con¬ 
ditions. 

As  a  first  attempt  to  deal  with  the  problem  of  tenuous  cirrus  clouds, 
collocated  IR  data  were  added  to  the  data  sample.  A  variable  threshold 
for  cirrus/no  cirrus  IR  count  was  selected  and  case  studies  begun  to  deter¬ 
mine  if  detection  of  cirrus  would  be  a  straightforward  procedure.  Analysis 
revealed  that  these  clouds  are  not  blackbody  radiators  and  varying  amounts 
of  background  radiation  leaks  through.  Therefore,  no  specific  IR  thresh- 
hold  can  be  considered  appropriate  for  even  a  majority  of  the  cases.  This 
threshold  value  must  be  interactively  determined  and  entered  into  the 
cloud  amount  algorithm  for  each  case.  High  cloud  studies  brought  to  light 
the  need  to  vary  the  size  of  the  analysis  box  with  varying  cloud  heights 
so  as  to  keep  the  satellite  field  of  view  in  agreement  with  that  of  a  sur¬ 
face-based  observer.  An  option  to  vary  the  size  of  the  analysis  box  has 
been  added  to  CENTER  but  remains  an  interactive  procedure. 

Threshold  techniques  for  specifying  cloud  cover  amounts  (Methods  3 
and  4)  are  a  more  accurate  method  than  average  brightness  techniques,  as 
is  shown  in  Table  1. 

By  properly  choosing  a  cloud/no  cloud  threshold,  with  a  decision- 
tree  or  similar  algorithm,  clear  and  overcast  conditions  can  be  classi¬ 
fied  to  a  high  degree  of  accuracy.  Sophistication  of  these  techniques 
can  help  to  enhance  the  results  achieved  by  the  simplistic  methods  test¬ 
ed  to  date. 


APPENDIX  A.  COMMAND  DOCUMENTATION  FOR  PROGRAM  CLUSTR. 


COMMAND:  CLUSTER 

PROGRAMMER:  G.  DENGEL 


KEY  IN :  CL 
DATE:  3  SEPT  81 

CHAPTER:  AFGL 


COMMAND  FORMAT: 


Cl.  OPT  PI  P2  P3  P4  P5  P6  P7 
COMMAND  DEFAULT  PARAMETERS:  CL 

cr 

Unless  an  option  is  specified,  CL  will  extract  the  VIS  and  IP  imagery 
samples  within  the  cursor,  cluster,  compute  histograms  of  texture  and 
gradient  for  each  cluster,  and  write  the  results  into  CLSFIL. 


PARAMETER  DESCRIPTION: 


OPT  Program  Options 


1 .  CL  LOAD  Area  Frame  Entry 

Area  =  digital  area  in  which  to  store  data 
Frame  =  TV  frame  on  which  to  display  data 
Entry  =  entry  number  in  CLSFIL  (^  10) 

2.  CL  TYPE  Entry  Cl  C2  C3  C4 

Assign  cloud  types  C1-C4  to  clusters  in  Entry 

3.  CL  HIS  Entry  Cluster  Function  Inc  Color  Max  X  Max  Y 


Draw  histogram  of  data  from  Entry  and  Cluster  specified. 

Function  =  TEX  or  GRD 
Inc  =1,  3,  or  5 

♦Color  =  1  (red) ,  2  (green) ,  3  (yellow) 

Negative  for  dashed  lines 
(default  =  green) 

Max  X  =  maximum  value  of  abscissa 
(default  =  100) 

Max  Y  -  maximum  value  of  ordinate 
(default  =  50) 


*  Axes  will  only  be  drawn  when  green  is  requested,  otherwise  an  overlay  is 
assumed. 


4.  CL  SUM  Ei  E2 . . .  (PRI) 

CL  SUM  ALL  (PRI) 

Print  data  summaries  for  entries 

El,  E2,  etc.  or  ALL  entries  in  file  on  CRT.  For  line  printer 
output,  last  parameter  should  be  "PRI." 

5.  CL  BOX 

Draws  boxes  around  all  samples  in  file  taken  from  currently 
displayed  frame. 
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6.  CL  SAVE  Tape  (NEW) 

Write  contents  of  file  to  magnetic  tape. 
Tape  =  Tape  number  (type  96) 

NEW  indicates  new  tape 

7.  CL  RES  Tape  File 

Restore  contents  of  tape  to  CLSFIL 
Tape  =  Tape  number  (type  96) 

File  =  Tape  file  to  be  restored 

8.  CL  INI 

nitialize  CLSFIL 

9.  CL  DEL  Entry 

Delete  entry  from  CLSFIL 


APPENDIX  B.  DOCUMENTATION  OF  THE  CLSFIL  FILE  STRUCTURE 

Type  Disk  File _ 

(Disk  File,  Savetape,  COMMON  block,  etc.) 

NAME  CLSFIL _ 

(If  applicable) 

DOCUMENTER  G.  Dengel _ _  Date  3  Sept  81 

Description  of  Structure: 

Type:  Data  file  (binary) 

Size:  1380  sectors 

Use:  Store  data  samples  and  results  for  program  CLUSTR. 

maximum  of  10  samples. 

Initialization:  Key  in  "CL  INI" 

Allocation:  138  sectors  per  entry 


Holds  a 


Imagery  samples  and  cluster  information: 


Sectors 


Words 


1-112 

VIS  data,  packed 

1-34 

VIS  data ,  packed 

35 

SSYYDDD 

36 

HHMMSS 

37 

Upper  left  corner 

line  coordinate 

38 

Upper  left  corner 

element  coordinate 

39 

Sample  number 

40 

Number  of  clusters 

41-43 

Cut  levels 

44-47 

Number  of  points  in  each  cluster 

48-51 

Mean  VIS  count  for 

each  cluster 

52-55 

Mean  IR  count  for 

each  cluster 

56-59 

Cloud  type  assigned  to  each  cluster 

60-81 

Cluster  frequency 

table,  packed 

82-111 

Unused 

112 

Number  of  entries 

in  file  (sector  16 

1-86 

IR  data,  packed 

87-112 

Unused 

Histogram  data,  half-words: 

Sectors 

Cluster 

Function 

18-22 

1 

TEX  1 

23-27 

1 

TEX  3 

28-32 

1 

TEX  5 

33-37 

1 

GRD  1 

38-42 

1 

GRD  3 

43-47 

1 

GRD  5 

48-77 

2 

78-107 

3 

108-137 

4 

A.  SUPPORT  TO  RESEARCH  PROJECTS 

The  McIDAS  Facility  supports  on  a  daily  basis  various  research  pro¬ 
jects  of  the  AFGL  Meteorology  Division.  The  Mesoscale  Forecast  Experi¬ 
ment  described  in  Section  2A,  for  example,  depends  wholly  on  the  McIDAS 
ability  to  display  meteorological  data.  Studies  in  satellite  meteorology 
are  particularly  dependent  on  McIDAS. 

In  addition,  McIDAS  is  used  to  support  field  programs.  During  1981 
eighteen  aircraft  flights  were  made  for  two  AFGL  research  projects:  the 
Large  Scale  Cloud  System  (LSCS)  and  Icing  studies.  Both  studies  used  a 
specially  instrumented  C-130  aircraft  based  at  Wright-Patterson  AFB,  OH. 

The  LSCS  study  sampled  large  homogeneous  cloud  masses  as  they  prog¬ 
ressed  from  west  to  east  across  the  U.S.  The  objective  was  to  obtain  a 
better  understanding  of  the  microphysics  of  clouds  in  large  winter  storms. 
The  hydrometeor  data  collected  by  the  aircraft  included  particle  size, 
number,  type,  distribution,  and  liquid  water  content.  The  Icing  study 
was  conducted  in  a  similar  manner  since  it  involved  observing  and  meas¬ 
uring  the  accumulation  of  ice  on  aircraft  control  surfaces.  The  C-130 
gathered  data  on  intensity  (light,  moderate,  heavy)  and  type  of  icing 
(glaze,  mixed,  rime) . 

McIDAS  satellite  imagery  loops  (both  visible  and  infrared)  depicting 
cloud  growth ?and  decay,  cloud  top  temperatures,  and  objective  analyses 
and  plots  of  radiosonde  data  (see  Fig.  1)  provided  the  ground  support  team 
information  on  the  best  sampling  locations  and  altitudes.  Additional  in¬ 
formation  examined  prior  to  and  in  real-time  during  the  flights  included: 
surface  observation  plots,  analyses  of  pressure  fields,  streamlines,  dew¬ 
point  temperatures,  wind  vectors,  precipitation  amounts,  current  weather, 
cloud  cover,  upper  air  cross-sections ,  icing  probability  (see  Table  1) , 
and  radar  observations  (see  Fig.  2) .  Data  for  both  studies  were  archived 
to  help  document  the  synoptic  situation  for  post  analysis. 

The  McIDAS  facility  was  also  used  to  collect  one-half  mile  resolu¬ 
tion  visible  satellite  imagery  last  winter  and  spring  over  North  America 
for  a  study  of  cloud  brightness  related  to  development  of  Tactical  Deci- 
sian  Aids.  Times  for  the  imagery  ranged  from  1530  to  1830  GMT  when  the 
sun's  declination  angle  was  at  its  highest  for  that  time  of  year.  Real¬ 
time  as  well  as  Sony  archived  imagery  were  used.  A  software  program  spe¬ 
cifically  written  for  the  project  provided  printouts  of  various  cloud 


Fig.  2.  Radar  Analysis  for  Eastern 


brightness  thresholds  found  in  the  imagery.  These  data  were  saved  on 
magnetic  tape  for  later  analysis. 


TABLE  1.  ICING  PROBABILITY  LISTING. 


*  *  *  ★ 

A  F  G 

L  M  C  I  D 

A  S 

*  *  *  * 

PROGRAM:  ICING 

ICING 

PROBABILITY  ESTIMATES  FOR 

72429  AT  10000 

Z  (FT) 

T  (0) 

TD  (0) 

T-TD 

TYPE 

INTENSITY  PROBABILITY 

1004 

.8 

-4.2 

5.0 

0 

*  1641 

0.0 

-4.8 

4.8 

0 

1880 

-.3 

-5.0 

4.7 

0 

1919 

0.0 

-5.3 

5.3 

0 

2389 

4.0 

-9.0 

13.0 

0 

3622 

3.8 

-10.2 

14.0 

0 

4731 

5.6 

-24.4 

30.0 

0 

5791 

5.6 

-24.4 

30.0 

0 

6890 

2.8 

-17.2 

20.0 

0 

7438 

2.6 

-7.4 

10.0 

0 

9289 

.  6 

-29.4 

30.0 

0 

*  9526 

0.0 

-22.4 

22.4 

0 

9882 

-.9 

-11.9 

11.0 

0 

10905 

-3.9 

-12.9 

9.0 

0 

11526 

-5.3 

-8.2 

2.9 

RIME 

LIGHT 

33 

11723 

-3.1 

-4.1 

1.0 

RIME 

MODERATE 

51 

13072 

-4.1 

-5.5 

1.4 

RIME 

MODERATE 

33 

14958 

-8. 1 

-9.0 

.  9 

MIXED 

MODERATE 

41 

18502 

-12.5 

-13.7 

1.2 

RIME 

MODERATE 

19 

23968 

-24.1 

-27.2 

3.1 

RIME 

LIGHT 

11 

*  INTERPOLATED  FREEZING  LEVEL 

McIDAS  provided  verification  data  for  a  project  by  the  Rome  Air 
Development  Center  (RADC/ET) .  This  study  examined  the  anomalous  prop¬ 
agation  of  microwave  siqnals  due  to  atmospheric  refractive  gradients. 

A  modified  McIDAS  software  program,  using  upper  air  data  from  the  Albany, 
NY  and  Chatham,  MA  radiosonde  stations,  provided  the  verification  data 
(Fig.  3  and  Table  2) . 
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ig.  3.  Refractive  Index  Plot. 


TABLE  2.  RADIOSONDE  DATA  LISTING. 


72317  011200 


984 

1.0 

.4 

070003 

281 

944 

-.4 

-.8 

115009 

609 

922 

-1.  3 

-1.6 

135013 

802 

908 

5.0 

4.7 

148015 

926 

876 

7.6 

7.1 

180023 

1219 

861 

9.0 

8.4 

188024 

1364 

850 

8.8 

8.0 

195025 

1471 

844 

8.6 

7.8 

190025 

1524 

814 

7.2 

6.5 

210025 

1829 

78- 

5.8 

5.2 

215024 

2134 

755 

4.3 

3.9 

215022 

2439 

727 

2.8 

2.  5 

220020 

2743 

700 

1.2 

1.0 

225019 

3058 

674 

-.9 

-10.4 

235023 

3353 

660 

-2.3 

-17.3 

237023 

3529 

649 

-2.6 

-13.6 

240023 

3658 

634 

-3.1 

-8.  1 

241022 

3848 

626 

-4.3 

-17.3 

241022 

3948 

612 

-4.5 

-11.5 

243021 

4126 

601 

-5.8 

-22.7 

245019 

4268 

587 

-7.5 

-37.5 

245019 

4453 

555 

-10.3 

-28.1 

245020 

4878 

539 

-11.9 

-22.9 

245020 

5111 

515 

-11.7 

-13.3 

245020 

5460 

500 

-11.9 

-13.5 

245019 

5687 

493 

-12.6 

-14.3 

245018 

5792 

473 

-14.6 

-16.8 

250017 

6097 

419 

-20.9 

-24.4 

275017 

7012 

400 

-23.3 

-27.4 

270020 

7357 

385 

-24.9 

-29.2 

260022 

7621 

369 

-26.9 

-31.2 

250026 

7926 

361 

-28.1 

-32.5 

251026 

8101 

311 

-36.4 

-42.1 

260031 

9146 

300 

-38.5 

-44.5 

260033 

9400 

250 

-48.9 

260048 

10624 

200 

-59.  3 

255051 

12054 

159 

-69.1 

263049 

13456 

150 

-70.5 

265048 

13803 

132 

-70.9 

270042 

14561 

127 

-66.9 

271040 

14792 

117 

-67.4 

275036 

15243 

100 

-68.7 

270030 

16229 

The  GOES  Research  Rapid  Scan  Winter 

Storm  project  was 

initiated  by 

AFGL 

,  in  November.  It 

calls  for 

five  winter  storm  research 

days  during 

the 

period  November  1981  -  January  1982. 

During  the  first 

storm  research 

day 

the  eastern  GOES 

satellite 

went  into 

a  seven-minute  scan  mode;  that 
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is,  it  transmitted  eight  images  in  one  hour  instead  of  the  normal  two. 
Areal  coverage  of  the  imagery  ranged  from  10°  N  to  80°  N  latitude, 
roughly  the  North  American  continent. 

The  rapid  scan  satellite  transmissions  were  recorded  on  the  MciDAS 
archive  system.  The  area  of  interest  was  then  stripped  off  and  stored  on 
magnetic  tape  for  future  display  and  analysis. 

B.  HARDWARE  ACQUISITION  AND  ENHANCEMENT 

The  process  for  storing  TV  images  in  MciDAS  used  to  display  satel¬ 
lite  imagery  was  changed  in  1981.  In  the  previous  system,  the  image  data 
flowed  from  the  host  computer  to  a  fast  storage  system  where  the  digital 
picture  was  converted  into  an  analog  signal  and  sent  to  the  analog  disk 
for  storage.  When  played  back  for  display,  the  picture  left  the  video 
disk  as  an  analog  signal,  entered  a  high  speed  analog  to  digital  converter 
where  it  was  converted  back  to  a  digital  image.  The  digital  data  were 
then  sent  to  the  red,  green  and  blue  enhancement  drawers.  After  enhance¬ 
ment  the  signal  was  converted  to  a  standard  RGB  analog  signal  and  sent  to 
the  video  monitors  for  display.  The  conversions  from  digital  to  analog 
and  back  again  contributed  to  picture  degradation. 

The  new  system,  installed  in  June,  is  the  Digital  Video  Storage  Sys¬ 
tem  (DVSS)  manufactured  by  Idetik  Corporation  of  Wisconsin.  This  system 
stores  the  image  as  sent  by  the  satellite,  in  digital  format,  thus  elim¬ 
inating  the  digital  to  analog  conversion.  The  system  is  composed  of:  a 
Harris  channel  interface  card,  two  V  ram  cards,  two  timing  and  control 
cards,  two  DVSS  cards,  two  86/12  processor  cards  with  on-board  memory,  two 
XY  logic  disk  controller  cards,  and  two  CDC  80  megabyte  Winchester  disks. 

The  DVSS  has  not  only  enhanced  picture  quality  but  has  also  increased 
the  picture  storage  volume  from  260  frames  per  terminal  to  400  frames  per 
terminal.  The  DVSS  was  installed  in  parallel  configuration  with  the  old 
analog  system  as  a  back-up.  Because  of  the  different  timing  outputs  of 
the  DVSS,  a  channel  interface  circuit  was  designed  to  accommodate  the  sig¬ 
nal  going  into  the  timing  drawer  for  Terminal  21.  Also  designed  into  this 
circuit  was  an  interlacing  capability;  that  is,  displaying  two  channels 
simultaneously. 

Another  component  purchased  from  and  installed  by  Idetik  Corporation 
was  a  VAS  r  ame  synchronizer.  This  unit  is  a  sophisticated  satellite  in- 


formation  processing  device  used  in  systems  receiving  Mode  A  and  Mode  AA 
transmission  formats.  It  is  designed  to  perform  major  and  minor  frame 
synchronization  on  stretched  visible  and  infrared  (VISSR  on  Mode  A)  and 
VISSR  atmospheric  sounder  (VAS,  or  Mode  AA)  data  from  the  new  GOES  satel¬ 
lites.  The  following  circuitry  cards  are  in  the  system:  Mode  A  frame 
synchronizer.  Mode  AA  frame  synchronizer,  hit  error  rate  and  line  count 
card,  GOES  satellite  data  simulator  card,  Intel  86/12  single  board  com¬ 
puter,  average  and  ingest  card,  arid  a  Multibus/Harris  channel  interface 
card.  When  the  McIDAS  data  ingest  and  handling  software  is  revised,  the 
frame  synchronizer  will  be  fully  operational. 

To  increase  McIDAS  operational  reliability  two  new  electronic  cir¬ 
cuits  were  designed  and  installed.  A  Texas  Instruments  Silent  700  elec¬ 
tronic  data  terminal  was  integrated  into  the  system  for  use  as  back-up 
for  the  McIDAS  teletype  terminal.  This  terminal  is  essential  for  commun¬ 
ication  between  the  McIDAS  operating  system  and  the  user.  The  interface 
circuit,  designed  by  Richard  Lynch  (AFGL) ,  can  accommodate  any  standard 
RS232  terminal  with  a  110  baud  rate.  The  other  circuit  installed  is  a 
relay  switching  panel  designed  to  reverse  the  controls  of  the  two  inter¬ 
active  consoles.  Using  coaxial  relays,  the  switching  panel  reverses  all 
RS232  lines  and  cables  to  and  from  the  joysticks  and  color  monitors. 

This  permits  what  is  displayed  on  one  terminal  to  be  viewed  on  the  other. 
Again,  this  is  a  back-up  measure  should  either  of  the  terminals  fail. 

During  the  year  the  elevation  motor  for  the  McIDAS  antenna  malfunc¬ 
tioned.  When  the  motor  was  sent  to  the  manufacturer's  maintenance  center 
for  repairs,  it  was  decided  to  overhaul  both  the  antenna  elevation  and 
azimuth  drive  assemblies.  A  new  motor  was  installed  in  the  azimuth  axis, 
all  drive  and  tachometer  belts  were  replaced,  and  all  gears  were  lubri¬ 
cated.  The  overhaul  was  performed  by  personnel  from  Scientif ic-Atlanta . 

The  Data  Products  line  printer  was  replaced  in  May  with  a  new  print¬ 
er  that  is  quieter  and  faster  (300  1pm  as  opposed  to  200  1pm)  than  its 
predecessor . 

C.  SOFTWARE  DOCUMENTATION 

Most  programs  written  for  the  AFGL  McIDAS  facility  rely  heavily  on 
pre-existing  subroutines  and  functions  that  are  on  disk  in  the  form  of 
library  modules.  A  main  library  and  seven  specialized  libraries  exist  in 


the  system.  The  main  library,  LL,  contains  a  wide  variety  of  support 
routines  for  functions  ranging  from  input/output  to  number  conversions. 

The  specialized  libraries  deal  with  the  areas  of  graphics,  nowcasting, 
navigation,  WINDCO,  real-time  radar,  manually  digitized  radar,  and  gen¬ 
eral  data  management.  Documentation  of  the  subprograms  within  these 
libraries  was  at  best  sparse,  consisting  of  occasional  comment  cards  in 
a  small  percentage  of  the  routines. 

SASC  completed  documentation  of  the  graphics  library  .,  GRALIB.  For 
each  of  the  58  modules  within  GRALIB,  the  type  of  subprogram,  calling 
sequence,  and  description  of  the  arguments  were  documented.  Any  addition¬ 
al  insights  into  use  of  the  routines  were  also  noted,  such  as  files  ac¬ 
cessed,  variable  types,  and  other  statements  that  may  be  helpful.  A  major 
index  of  all  graphics  library  modules  was  also  compiled,  giving  the  name 
of  each  routine  and  a  short  description  of  its  purpose. 


SECTION  5. 
RADAR  METEOROLOGY 
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A.  WIND  RETRIEVAL 


I.  Introduction 

Several  methods  for  retrieving  winds  from  single-Doppler  radar  have 

been  proposed.  Browning  and  Wexler  (1968)^"  and  Waldteufel  and  Corbin 
2 

(1979)  presented  several  such  techniques  which  depend  upon  linearity 
of  the  wind  velocity  components  and  resolve  the  wind  over  large  areas. 

Various  storm  and  echo  tracking  schemes  have  also  been  presented  (e.j., 

3  4  5 

Crane,  1979  ,  and  Bjerkaas  and  Forsyth,  1979  ).  Rinehart  (1979)  and 

Smythe  (1981) ^  presented  techniques  which  give  local  estimates  of  t;, 
wind  over  much  stiller  areas  by  seeking  smaller-scale  patterns  in  reflec¬ 
tivity  fields  usinq  correlation.  Smythe  (1981)  also  investigated  the 
feasibility  of  tracking  small-scale  perturbations  in  Doppler  radial  ve¬ 
locity  fields  to  estimate  the  horizontal  wind  and  produced  favorable  re¬ 
sults  when  these  winds  were  compared  with  winds  synthesized  from  two 
Doppler  radars.  This  report  is  a  continuation  of  the  study  of  Smythe 
(1981) . 

II.  Discussion 

The  algorithm  of  Smythe  (1981)  was  modified  to  increase  its  effi¬ 
ciency  in  an  attempt  to  meet  AFGL  requirements  for  real-time  applications. 
Correlation  coefficients  were  computed  using  the  equation 
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where  r  is  the  coefficient  of  cross-correlation,  the  x. 's  and  v.'s  are 
either  reflectivity  or  pulse  pair  velocity  estimates  for  each  pulse  vol¬ 
ume  in  the  regions  correlated,  and  N  is  the  number  of  samples  available, 
as  some  data  might  have  been  removed  due  to  various  editing  criteria. 

Doppler  velocities  were  averaged  over  the  correlation  region  (array)  of 
the  first  scan.  This  average  was  used  to  find  a  first-guess  approxima¬ 
tion  of  the  radial  displacement  of  the  region  at  a  later  time.  Data  (re¬ 
flectivity  or  velocity)  were  then  lagged  azimuthally  to  search  for  a  cor¬ 
relation  maximum.  Fig.  1  shows  an  example  of  the  first-guess  displace¬ 
ment  (D) ,  the  array  of  the  first  scan  (hatched  area) ,  and  the  region  of 
the  second  scan  (between  dashed  curves)  over  which  the  initial  array  was 
lagged  and  correlated. 

Several  software  versions  of  the  wind  retrieval  algorithm  were  writ¬ 
ten,  two  of  which  were  hybrids  consisting  of  both  Assembler  and  Fortran 
statements.  Portions  of  the  programs  were  timed  and  the  tightest  loops 
were  recoded  from  Fortran  to  Assembler.  In  these  loops  floating  point 
arithmetic  was  changed  to  halfword  arithmetic.  Also,  register  operations 
were  used  extensively.  Excessive  iterations  through  unnecessary  time- 
consuming  code  were  thus  avoided. 

Using  the  hybrids  it  was  found  that  the  time  required  to  analyze  an 
area  of  dimensions  40  km  x  40°  (which  is  large  enough  to  cover  a  thunder¬ 
storm)  was  approximately  140  s  on  the  Interdata  7/32  minicomputer.  Although 
the  software  was  not  tested  on  the  Perkin-Elmer  (P-E)  3242,  run-time  should 
be  less  than  one-third  of  the  present  run-time.  This  would  then  allow  a 
real-time  analysis  of  six  levels  within  the  volume  of  the  storm. 

Features  not  coded  in  the  original  programs  were  added.  Those  features 
included  the  flagging,  before  processing,  of  arrays  containing  high  per¬ 
centages  of  low  signal  data,  thus  distinguishing  areas  of  storm  (strong 
echo  return)  from  weak  signal.  The  code  was  also  modified  to  accommodate 
reflectivity  and  velocity  data  staggered  at  different  range  gate  intervals 
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Fig.  1.  Wind  Retrieval  Technique  Geometry 


As  (e.g.,  As  .  150m,  As  .  .  . ^  =  300m). 

velocity  reflectivity 

Shown  in  Fig.  2  are  correlation  coefficients  computed  from  data  col- 

o 

lected  May  2,  1978,  at  0.6  elevation  m  Oklahoma.  The  coefficients  were 
obtained  from  correlating  reflectivity  and  lagging  in  intervals  of  0.6° 
azimuthally  and  300m  radially  about  the  first-guess  approximation  of  the 
radial  displacement.  The  temporal  lag  was  about  min.  Comparable 
maxima  m  orrelation  coefficients  are  at  lags  -8  and  +6.  In  this 
example  ine  structure  of  the  field  appears  to  be  wave-like.  Fig.  3  il¬ 
lustrates  this  wave  structure  also.  Smythe  (1981)  showed  that  this  type 
of  structure  can  obscure  the  true  maximum  and  produce  spurious  results. 
Further  investigation  is  required  to  identify  these  features  and  deter¬ 
mine  causal  mechanisms. 

Fig.  4  depicts  a  wind  field  reconstructed  from  May  2,  1978  data  col¬ 
lected  at  0.7°  elevation,  for  which  3  km  x  4°  correlation  arrays  were 
used.  Light  ends  of  arrows  correspond  to  arrowheads.  Shown  is  a  gust 
front  marking  the  leading  edge  of  cold  air  outflow  ahead  of  a  storm  ad¬ 
vancing  toward  the  radar.  Convergence  is  clearly  evident. 

An  algorithm  for  reducing  effects  of  range  ambiguities  in  velocity 
data  was  also  added.  If  the  pulse  width  of  reflectivity  were  expanded 
to  double  the  width  of  velocity  sampling,  then  velocity  values  can  be  as¬ 
signed  to  their  respective  trips.  Reflectivity  (power)  values  at  range 
R  were  compared  with  values  at  range  2R,  If  the  difference  in  power  was 
greater  than  7  dB  the  range  with  the  larger  power  was  assigned  the  velo¬ 
city  value,  the  other  was  flagged.  If  the  difference  was  less  than  7  dB 
neither  range  was  assigned  a  velocity  value. 

Another  approach  to  the  problem  is  to  use  two  elevations  in  the  anal¬ 
ysis.  It  was  thought  that  by  this  approach  a  data  volume  might  be  recog¬ 
nized  more  readily  than  if  only  one  level  were  used.  This  could  also 
lead  to  reconstructing  three-dimensional  winds.  A  program  to  do  this  an¬ 
alysis  was  coded  but  not  tested. 

At  the  end  of  the  contract  period  cases  of  April  28,  1978  and  April  10, 
1979  were  being  examined.  In  the  1978  case  data  were  collected  in  volume 
scan  mode  with  full  plan  position  indicators  (PPI's)  ,  so  scans  at  the  same 
level  were  a±>out  six  minutes  apart.  However,  in  the  1979  case  data  were 
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Fia.  2.  Correlations  fron  Scans  Collected  Gi  Ilinutes  Apart  of  an  Area  of 
Storm  on  tiay  2,  197G. 


collected  at  0.5°,  so  scans  were  approximately  one  minute  apart.  The  1979 
data,  collected  by  the  National  Severe  Storms  Laboratory  (NSSL)  but  ob¬ 
tained  through  the  courtesy  of  the  National  Center  for  Atmospheric  Re¬ 
search,  include  dual-Doppler  data.  For  the  other  cases  no  dual  radar  data 
were  available  for  comparisons. 

Since  the  NSSL  data  were  packed  into  six-bit  words,  it  was  necessary 
to  write  software  to  upack  these  data  so  that  they  would  conform  to  the 
32-bit  word  format  of  the  AFGL  P-E  computers.  This  software  was  combined 
with  editing  software. 

For  analyses  of  the  April  28,  1978  and  April  10,  1979  data  the  algor¬ 
ithm  of  Smythe  (1981)  was  changed  so  the  the  range  of  uncertainty  about 
the  first-guess  range  displacement  was  not  fixed  for  each  correlation  ar¬ 
ray  but  instead  was  allowed  to  vary  according  to  the  standard  deviation 
of  level  of  confidence  limits  about  the  mean  velocity.  Preliminary  re¬ 
sults  of  testing  were  encouraging. 

B.  REAL-TIME  STORM  CELL  DETECTION  AND  TRACKING 

SASC  continued  its  efforts  to  develop  a  real-time  software  technique 
for  the  detection  and  tracking  of  significant  storm  cells  by  reducing  the 
running  time  of  the  processing  algorithm  originally  developed  by  Crane 
(1979) .  It  extracts  the  essential  information  required  to  describe  the 
organization,  development,  and  motion  of  the  storm  system.  This  informa¬ 
tion  includes  the  characteristics  of  the  active  cells  within  the  storm 
system,  the  clusters  of  cells,  the  f ixed-threshold  contours  or  echo  re¬ 
gions  that  encompass  the  cells,  and  the  complexes  of  echo  regions  that 
comprise  the  storm  system.  The  algorithm  was  modified  in  1979  for  use  on 
the  Interdata  7/32  on-site  computer  at  AFGL. 

The  original  software  includes  one  main  program,  nine  subroutines, 
and  one  data  block.  The  main  program  and  the  data  acquisition  subrou¬ 
tine  are  written  in  a  mixed  mode  of  Fortran  and  assembler;  other  subrou¬ 
tines  are  written  in  Fortran.  The  software  is  designed  to  collect  raw 
weather  radar  data,  decode,  perform  the  necessary  data  extraction,  and 
provide  these  data  for  later  analysis.  To  make  this  system  real-time 
the  data  must  be  gathered,  processed,  and  stored  quickly  enough  to  re¬ 
ceive  the  next  set  of  input. 

The  data  processing  scheme  is  shown  in  Fig.  5.  The  data  decoding 
and  conversion  processing  routine  in  the  data  acquisition  subroutine  is 


Fig.  5.  Overall  Processing  Scheme. 
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the  most  time-critical  process.  For  instance,  if  there  are  768  range 
cells  per  radial  and  the  azimuthal  difference  (A0)  is  0.8  degree,  this 
routine  needs  to  be  iterated  more  than  300,000  times.  Event  ilcula- 
tion  and  referenced  peak  calculation  routines  are  next  most  pressing 
in  the  time  constraint  problem.  They  have  to  be  run  over  400  times  per 
azimuth  scan  if  the  azimuthal  difference  is  0.8  degree. 

Since  assembly  language  is  capable  of  using  registers  directly  for 
central  processor  unit  (CPU)  execution,  it  is  more  efficient  in  bit  man¬ 
ipulation  and  addressing  calculation,  which  are  the  important  steps  in 
reducing  CPU  execution  time.  The  conversion  and  decoding  processing 
routine  and  event  calculation  routine  were  therefore  rewritten  into  as¬ 
sembly  language.  In  addition,  the  conversion  algorithm  from  bit-count  to 
dBZ  was  optimized.  These  modifications  saved  more  than  50  sec  per  azi¬ 
muth  scan. 

On  the  other  hand,  assembler  is  hard  to  maintain  and  needs  more  de¬ 
bugging  time.  Also,  a  mixed  mode  of  Fortran  and  assembler  needs  extra 
housekeeping  time  for  interchange  between  the  two  different  languages. 
Therefore  the  Fortran  code  was  optimized  in  the  referenced  peak  calcu¬ 
lation  subroutine. 

April  30,  1978  data  were  employed  to  test  these  modifications.  The 
output  of  the  post  mission  (1979)  program  was  used  as  the  standard. 

Timing  analysis  of  the  modified  version  of  the  tracking  program  showed 
that  the  average  processing  time  for  each  azimuth  scan  on  the  Interdata 
7/32  is  about  47  sec,  which  is  close  to  the  requirement  for  real-time 
operation  (40  sec) . 

However,  a  problem  was  detected  with  the  post  mission  output:  the 
same  program,  using  identical  input  data,  on  successive  runs  on  the  Inter 
data  7/32  minicomputer  produced  two  different  outputs.  Results  for  one 
run  are  shown  in  Fig.  6,  which  is  a  plan  view  of  the  centroids  of  the 
cells  detected  by  the  software.  On  a  subsequent  run  the  four  cells  ident 
ified  by  arrows  were  not  detected  at  those  locations.  Instead,  four  cell 
were  detected  at  large  distances  north  or  south  (off-scale  on  the  plot) , 
where  none  had  been  previously  detected.  The  source  of  this  problem  has 
not  been  identified.  As  a  consequence,  use  of  the  post  mission  program 
as  a  standard  must  be  questioned.  There  is,  however,  agreement  among  the 
the  results  shown  in  Fig.  6,  the  results  from  an  earlier  run  made  with 
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the  same  program,  and  the  results  from  the  current  modified  program.  This 
consistency  suggests  that  the  changes  incorporated  in  the  modified  ver¬ 
sion  are  valid  but  that  the  two  programs  may  share  a  common  problem. 

C.  TURBULENCE  MEASUREMENTS 

Several  tasks  were  undertaken  in  support  of  AFGL's  Wallops  Island 
turbulence  project.  A  program  was  written  to  read  data  from  fourteen- 
track  tape  via  the  pulse-pair  recorder  interface  and  write  the  data  to 
nine-track  tape.  Several  nine-track  data  tapes  were  prepared  for  future 
use  in  non-real-time  analyses.  Coding  of  software  for  unpacking  raw  data 
into  16-bit  halfwords  was  begun.  Also,  writing  of  a  program  tc  produce 
B-scans  of  data  was  begun. 

A  study  based  on  techniques  proposed  by  Atlas  and  Srivastava  (1969)  ' 

8 

and  Serafin  (1972)  was  initiated.  Implementation  of  these  techniques 
would  permit  estimation  of  eddy  dissipation  rates  by  incoherent  (non-Dop¬ 
pler)  radar. 
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